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We discuss an extension of the fluctuation theorem to stochastic models that, in the limit of zero 
external drive, are not able to equilibrate with their environment, extending results presented by 
Sellitto (cond-mat/9809186). We show that if the entropy production rate is suitably defined, its 
probability distribution function verifies the Fluctuation Relation with the ambient temperature re- 
placed by a (frequency-dependent) effective temperature. We derive modified Green-Kubo relations. 
We illustrate these results with the simple example of an oscillator coupled to a nonequilibrium bath 
driven by an external force. We discuss the relevance of our results for driven glasses and the diffu- 
sion of Brownian particles in out of equilibrium media and propose a concrete experimental strategy 
to measure the low frequency value of the effective temperature using the fluctuations of the work 
done by an ac conservative field. We compare our results to related ones that appeared in the 
literature recently. 



I. INTRODUCTION 



Relatively few generic results for non-equilibrium systems exist. Recently, two such results that apply to seemingly 
very different physical situations have been proposed and extensively studied. One is the fluctuation theorem that 
characterizes the fluctuations of the entropy production over long time-intervals in certain driven steady states [1, 2]. 
Another one is the extension of the fluctuation-dissipation theorem that relates induced and spontaneous fluctuations 
in equilibrium to the non-equilibrium slow relaxation of glassy systems [3, 4]. While the former result has been proven 
for reversible hyperbolic dynamical systems [2, 5] and for the driven stochastic dynamic evolution of an open system 
coupled to an external environment [6, 7], the latter has only been obtained in a number of solvable mean-field models 
and numerically in some more realistic glassy systems (see [8] for a review). The modification of the fluctuation- 
dissipation theorem can be rationalized in terms of the generation of an effective temperature [9, 10]. The expected 
thermodynamic properties of the effective temperature have been demonstrated in a number of cases [11]. 

One may naturally wonder whether these two quite generic results may be included in a common, more general 
statement. The scope of this article is to discuss this possibility in general, illustrating it with a very simple example 
in which one can very easily reach the 'driven limit' and the 'non-equilibrium relaxational' case. This project was 
pioneered by Sellitto [12] who asked the same question some years ago and tried to give it an answer using a stochastic 
lattice gas with reversible kinetic constraints in diffusive contact with two particle reservoirs at different chemical 
potentials. Other developments in similar directions have been proposed and analyzed by several authors [13-15]. We 
shall discuss them in Sect. VIII. 

Before entering into the details of our calculations let us start by reviewing the precise statement of the usual 
fluctuation theorem and the extension of the fluctuation-dissipation theorem as well as the definition of the effective 
temperature. 
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A. The fluctuation theorem 



The fluctuation theorem concerns the fluctuations of the entropy production rate [16], that we call a(x), in the 
stationary non-equilibrium state of a dynamical system defined by a state variable x <G M, where M is the phase 
space of the system. 

In non-equilibrium stationary states the function a(x) has a positive average, a+ = J M dxfj, + (x)a(x) > 0, where 
yu + (a;) is the stationary state distribution. This allows one to define 

i i r' 2 

p(x) = S T = / dt o-(S t x) , (1) 

TCT+ TCT+ J_ r/2 

where S t x, with time t G [— t/2,t/2], is a segment of the system's trajectory starting at the point x at t = — r/2. 
The large deviation function of p(x) is 

C(p) = lim r-HniTrip) , (2) 

where n T (p) is the distribution of p in the stationary non-equilibrium state. The fluctuation theorem is the following 
statement about the large deviation function of p(x): 

c(p) - a-p) = po- + . (3) 

The relation (3) was first discovered in a numerical simulation [1] , and subsequently stated as a theorem for reversible 
hyperbolic dynamical systems (the Gallavotti- Cohen Fluctuation Theorem) [2, 5]. The proof was extended to Langevin 
systems coupled to a white bath in [6] and to generic Markov processes in [7]. Equation (3) was successfully tested 
in a wide number of numerical simulations, see e.g. [17-20], and, more recently, it was also analyzed in experiments 
on granular materials and turbulent flows [21-24]. At present, it is believed to yield a very general characterization 
of the fluctuations of the entropy production in a variety of out of equilibrium stationary states. 

In all the cases cited above an external drive maintains the systems in a stationary non-equilibrium regime. In the 
absence of the drive the systems so far analyzed easily equilibrate. In the stochastic case, the systems are in contact 
with an equilibrated environment at a well-defined temperature. 



B. The effective temperature 

The analytic solution to the relaxation of mean-field glassy models following a quench into their glassy phase 
demonstrated that their relaxation occurs out of equilibrium [3, 4]. The reason why these models do not reach 
equilibrium when relaxing from a random initial condition is that their equilibration time diverges with their size. 
Thus, when the thermodynamic limit is taken at the outset of the calculation, all times considered are finite with 
respect to the size of the system. These systems approach a slow nonequilibrium regime in which one observes a 
breakdown of stationarity and, more importantly for the subject of this paper, a violation of the fluctuation-dissipation 
theorem that relates spontaneous and induced fluctuations in thermal equilibrium. 

The relation between spontaneous and induced fluctuations found in mean-field glassy models is, however, sur- 
prisingly simple. Let us define the linear response of a generic observable O measured at time t to an infinitesimal 
perturbation constantly applied since a previous 'waiting-time' t w , and the correlation between the same (unperturbed) 
observable measured at t,„ and i, 



t r t 

(4) 

h=0 



X (t,t w ) = I dt' R(t,t') = I dt' S (°^ 



5h(f) 

C(t,t w ) = (0(t)0(t w )) , (5) 



where, for simplicity, we assumed that the observable O has a vanishing average, (0(t)) = for all t. In the cases we 
shall be interested in the relation between these two quantities takes the form 

lim X (t,t w ) = f[C(t,t w )} (6) 

in the long waiting-time limit after the initial time to- This equation means that the waiting-time and total time 
dependence in x enters only through the value of the associated correlation between these times. This is trivially true 
in equilibrium since the fluctuation- dissipation theorem states 

/(C) = i(l-C), (7) 
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for all times t > t w >t eq , where T is the temperature of the thermal bath (and the one of the system as well) and we 
set the Boltzmann constant to one, fc# = 1. Out of equilibrium /(C) can take a different form. In a large variety of 
models with slow dynamics that describe aspects of glasses /(C) is a broken line, 

/(C) = i(l-C) 6(C -q ea ) + 

see [8] for a review. This broken line has two slopes, — 1/T for C > q ea (i.e., small t — t w ), and —1/T e ff for C < q ea 
(i.e., large t — t w ). The breaking-point q ea has an interpretation in terms of intra-cage and out-of-the cage motion in 
the relaxation of structural glasses, inter-domain and domain wall motion in coarsening systems, etc. Since T e ff is 
found to be larger than T the second term violates the fluctuation-dissipation theorem. We have used the suggestive 
name effective temperature, T e jj, to parametrize the second slope. The justification is that for mean- field glassy 
models — and within all resummation schemes applied to realistic ones as well — T e ff does indeed behave as a 
temperature [9] . A similar relation was later found numerically in the slow relaxational dynamics of a number of more 
realistic glassy systems such as Lennard- Jones mixtures [25]. More complicated forms, with a sequence of segments 
with different slopes appear in mean-field glassy models of the Sherrington-Kirkpatrick type. 
The general definition of the effective temperature is 

Teff(C) do' K> 

In order to be consistent with the thermodynamic properties one needs to find a single value of T e ff in each time 
regime as defined by the correlation scales of Ref. [4] . 

So far we considered systems relaxing out of equilibrium in the absence of external drive. The extension of the 
definition of the effective temperature to the driven dynamics of glasses and super-cooled liquids was again motivated 
by the analytic solution of mean- field models under non-potential forces [26-28] . The existence of non-trivial effective 
temperatures, i.e. of functions /(C) that differ from Eq. (7), was also observed in numerical simulations of sheared 
Lennard-Jones mixtures [29, 30] and in a number of other driven low-dimensional models [31]. 

In the case of relaxing glasses the dynamics occurs out of equilibrium because below some temperature the equi- 
libration time falls beyond all experimentally accessible time-scales. These macroscopic systems then evolve out of 
equilibrium even if they are in contact with a thermal reservoir, itself in equilibrium at a given temperature T. Un- 
der the effect of stirring forces, supercooled liquids and glasses are typically driven into a nonequilibrium stationary 
regime, even for relatively modest flow. 

C. A connection between the two? 

The fluctuation theorem and the fluctuation-dissipation theorem are related: indeed, it was proven that, for 
systems which are able to equilibrate in the small entropy production limit (er + — > 0), the fluctuation theorem 
implies the Green-Kubo relations for transport coefficients, that are a particular instance of the fluctuation-dissipation 
theorem [6, 32, 33]. 

It is now natural to wonder what is the fate of the fluctuation theorem if the system on which the driving force 
is applied cannot equilibrate with its environment and evolves out of equilibrium even in the absence of the external 
drive. Our aim is to investigate whether the fluctuation theorem is modified and, more precisely, whether the effective 
temperature [9] enters its modified version. In particular, this question will arise if the limit of large sampling time, 
t in Eqs. (1) and (2), is taken after the limit of large system size. The order of the limits is important because a 
finite size undriven system will always equilibrate with the thermal bath in a large enough time r. As the fluctuation 
theorem concerns the fluctuations of a for r — > oo, if one wants to observe intrinsic nonequilibrium effects, the latter 
limit has to be taken after the thermodynamic limit. Some conjectures about this issue have recently appeared in the 
literature [12-15] and we discuss them in Sect. VIII. 

Our idea is to study the relaxational and driven dynamics of simple systems such that the effective temperature is 
not trivially equal to the ambient temperature. For a system coupled to a single thermal bath, this happens whenever: 

(i) the thermal bath has temperature T, but the system is not able to equilibrate with the bath. This is realized 
by the glassy cases discussed above, provided the sampling time is smaller than the equilibration time; and/or 

(ii) the system is very simple (not glassy) but it is set in contact with a bath that is not in equilibrium. One can 
think of two ways of realizing this. One is with a single bath represented by a thermal noise and a memory friction 
kernel that do not verify the fluctuation-dissipation relation [34]. This situation arises if one considers the diffusion 
of a Brownian particle in a complex medium (e.g. a glass, or granular matter) [35-37]. In this case the medium, 
which acts as a thermal bath with respect to the Brownian particle, is itself out of equilibrium. Another example 



id 



Qea) + Tf, (lea ~ C) 



L eff 



0(q ea - C) 



(8) 
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is the one of a system coupled to a number of equilibrated thermal baths with different time-scales and at different 
temperatures [10]. 



The cases (i) and (it) mentioned in the previous section are closely related as, at least at the mean-field level, the 
problem of glassy dynamics can be mapped onto the problem of a single "effective" degree of freedom moving in an 
out of equilibrium self-consistent environment [10, 34]. Situations (i) and (it) are then described by the same kind of 
equation, namely, a Langevin equation for a single degree of freedom coupled to a non-equilibrium bath. 

Indeed, in the study of mean-field models for glassy dynamics [10, 34] and when using resummation techniques within 
a perturbative approach to microscopic glassy models with no disorder, it is possible to reduce the TV-dimensional 
equations of motion to a single equation for an 'effective' variable by means of a saddle-point evaluation of the dynamic 
generating functional. These equations are valid in the large size N — > oo limit; the discussion that follows concerns 
the fluctuation relation involving S T for r — > oo taken after N — > oo: for these finite timescales r a fluctuation relation 
involving effective temperatures may and will arise, while for the extreme times r — ► oo before N — > oo the usual 
fluctuation theorem (involving only the bath temperature T) holds. 

The effective equation of motion of a single spin at finite time reads 



where pt is a Gaussian noise such that (ptPf) — D(t, t'). The "self-energy" E(i, t') and the "vertex" D(t, t') depend on 
the interactions in the system. They cannot be calculated exactly for generic interactions but they can be approximated 
within different resummation schemes (mode-coupling, self-consistent screening, etc.) or calculated explicitly for 
disordered mean-field models. The last term (if present) represents an external drive. 

In the absence of the driving force Eq. (10) shows a dynamical transition at a temperature T<j. Above Td, in the 
limit t, t' — ► oo the functions E(i, t') and D(t, t') become time-translation invariant and are related by the fluctuation- 
dissipation relation with E playing the role of a response and D being a correlation. Below Td the system is no more 
able to equilibrate with the thermal bath and ages indefinitely, i.e. the functions E(i, t') and D(t, t') depend on t and 
t' separately even in the infinite time limit, and the relaxation time r a grows indefinitely. At low temperatures the 
fluctuation-dissipation theorem does not necessarily hold, and the relation between E and D can be used to measure 
the effective temperature for the particular problem at hand. 

In the case of a driven mean field system [26-28], the external force is also present in Eq. (10) and after a transient 
the system becomes stationary for any temperature, i.e. p{t) = p, T,(t,t') = E(t — t'), and D(t,t') — D(t — t'). The 
functions D and E depend on the strength a of the driving force. Below Td and for small a they are again related by 
a generalization of the fluctuation-dissipation theorem similar to the one obtained for a = 0. 

In this paper we focus on Langevin equations similar to Eq. (10). Their main characteristic is that the thermal 
bath, represented by the functions T,(t,t') and D(t,t'), is not in equilibrium. As we are interested in the driven case, 
we restrict to considering stationary functions E(i — t') and D(t — t'). 

Let us remark that Eq. (10) is expected to describe the dynamics of a single Brownian particle immersed in a 
non-equilibrium environment, or the dynamics of an effective degree of freedom representing a many-particle mean- 
field glassy system. If one wishes to describe in full detail the relaxation of real glasses in finite dimensions, more 
complicated effects have to be taken into account. For example, the dynamics is expected to be heterogeneous yielding 
a local effective temperature which may depend on space inside the sample, sec e.g. [38] for a detailed discussion. The 
extension of the results that we shall present in this paper to glassy systems in finite dimension might require additional 
work. 



The aim of this work is to discuss the validity of the fluctuation theorem for the Langevin equation (10) in presence 
of a non-trivial environment represented by the functions E(£) and D(t). This equation describes a situation in which 
the relaxing system does not equilibrate with its environment (even in the absence of driving forces) and a non-trivial 
effective temperature defined from the modification of the fluctuation-dissipation theorem exists. Our first aim is to 
identify the entropy production rate and to show that the effective temperature of the environment should replace 
the ambient temperature of a conventional bath. 

As particular cases we investigate analytically the dynamics in a harmonic well and numerically a case in which the 
equations of motion are nonlinear. The former problem is relevant for experiments on confined Brownian particles in 



D. Effective equations for the dynamics of mean- field glasses 




(10) 



E. Summary and structure of the paper 
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complex media [35, 39]. In both cases we show that we show that the entropy production rate that we introduced 
verifies the fluctuation relation in general. 

The paper is organized as follows. In section II we identify the correct definition of entropy production rate using a 
rather general procedure proposed by Lebowitz and Spohn. As an illustration, in section III we analytically compute 
the large deviation function in the harmonic case for a white equilibrium bath, already discussed in Ref. [6], and for 
a complex bath. We show that the large deviation function is a convex function of p and satisfies the fluctuation 
relation. In section IV we compare the analytic results obtained for the linear problem with numerical data. We also 
numerically investigate a nonlinear Langevin equation, in which the entropy production is not only given by the work 
of the external forces, but should contain an 'internal' term as well. Surprisingly, this term turns out to be negligible, 
a result we attribute to decorrelation between the work of internal and external forces. In section V we make contact 
with glassy problems and discuss some connections with recent numerical simulations. Section VI is devoted to the 
analysis of the link between the modified fluctuation theorem and modified Grccn-Kubo relations. In Section VII we 
show cxplicitely that the distribution of the work done by a slow periodic drive satisfies the fluctuation relation with 
the low-frequency effective temperature, and argue that this procedure is the one that could be "easily" implemented 
experimentally as a means to measure the (low frequency) effective tempeature. In the conclusions we briefly discuss 
some experiments that could test our predictions, and compare them to previous studies of the same problem. 



II. ENTROPY PRODUCTION 



In this Section we introduce the set of dynamic models that we discuss in detail, namely, stochastic processes of 
generic Langevin type with additive noise. Next, we define the large deviation function that we shall use to test the 
validity of the fluctuation relation and, finally, we derive a general expression for the entropy production rate. 



A. The model 



In this article we focus on different aspects of the random motion of a particle in a confining potential, in contact 
with a thermal environment, and under the effect of a driving external force. The Langevin equation describing the 
motion of such a particle in a d dimensional space reads 

f°° dV 
mr a (t) + dt'g afj (t-t')r p (t') = - — [r(t)]+p a (t) + h a (t), a=l,...,d. (11) 

J — oo or a 

r = (r\, . . . , r^) is the position of the particle. We pay special attention to the case d ~ 2 and call (x,y) the two 
components of the position vector, m is the mass of the particle and V(r) is a potential energy. As an example 
we shall work out in detail the simple harmonic case, V(r) = -| ^2 a r 2 a with k the spring constant of the quadratic 
potential. p(t) is a Gaussian thermal noise with zero average and generic stationary correlation 

(Pa{t)p {t'))=5 a pv(t-t') a,0 = l,...,d, (12) 

with v{t — t') a symmetric function of t — t' . The memory kernel g a fi{t — t') extends the notion of friction to a more 
generic case. We assume a simple spatial structure, g a p(t — t') = S a fj g(t — t 1 ). In order to ensure causality we take 
g(t — t') to be proportional to 9(t — t'). We define 

g(t) = 6(t)f(t) with f(t) = g(t)+g(-t) . (13) 

The initial time to has been taken to — oo. h{t) is a time-dependent field that we either use to compute the linear 
response or represents the external forcing. 

It will be useful to introduce Fourier transforms. We use the conventions 

/OO J poo 
e--VM, pM=/ dte^p(t). (14) 

-oo ^ J -oo 

The fluctuation-dissipation theorem states that for systems evolving in thermal equilibrium with their equilibrated 
environment the linear response is related to the correlation function of the same observable as 

m -m.*®., sfl =m - t) -«(*,], da. 
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see Eqs. (4) and (5) for the definitions of R and C. Here and in what follows we set the Boltzmann constant to one, 
k B = 1. To write the second expression we used that C(t) is an even function of t (and C(t) is odd) and defined 
9(0) = 1/2 (the same convention is used in the rest of the paper). After Fourier transforming the second expression 
becomes 

loC(lo) = 2TlmR(uj) , (16) 

and the real part of R(u>) is related to ImR(uj) by the Kramers-Kronig relation. 

The functions g and v are the integrated response and correlation of the bath, respectively. As will be discussed in 
detail in section III, if the bath is itself in equilibrium at a temperature T they are related as 



1 _ 2Reg(oj) 
T ~ v(u) ' 

For an out of equilibrium bath we define the frequency dependent temperature 

1 2Reg(uj) 



T(lu) ' v(uj) 



(17) 



(18) 



and its inverse Fourier transform 



T-\t) = I ^L^9M e -i^t = I dt > u -i {t tV{t > } . (19) 



2tt y(u) 



where 



v -\ t) = / ^_L e -*-* . (20) 

w J 2ir u(u) V ' 

Note that both f(t) and T _1 (t) are even functions of t. If the bath is in equilibrium at temperature T, T _1 (t — t') = 
S(t — t')/T. We shall assume throughout that T~ 1 (t) goes to zero fast enough for large t. 

The main result of this paper is that the fluctuation relation for the probability distribution function tt t (p) holds 
for long times r: 

ra + p ~ In (21) 



with 

r/2 r r/2 



S T = T(j + p= dt I dt'T-^t-t') r a (t) 

J-t/2 J-t/2 



(22) 



B. Large deviation function 



In the rest of this section we identify the entropy production rate for Eq. (11) following the procedure proposed by 
Lebowitz and Spohn [7]. 

The fluctuation relation is a symmetry property of the probability distribution function (pdf) of the entropy 
production rate at that can be also expressed as a symmetry of the large deviation function. Calling a + the average 
value of at, the PDF of the variable p defined in Eq. (1) is defined as 



TT T ( P )=V 



ra+ 



= P 



= Ta + (8(S T -pra + )) 



(23) 



where the angular brackets denote an average over the realizations of the noise. The large deviation function [nor- 
malized so that C(l) = at the maximum] is given by 



C(p) = lim r- 1 ln[7r T (p)/7r T (l)] . 

T— >00 

It is easier to compute the characteristic function 

0(A) = — lim r _1 ln(cxp[— \a T \) ; 



(24) 



(25) 
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the latter being the Legendre transform of ((p). Indeed, 

f r 7„ p T[C(p)-Ap ( T + ] T max p [C(p)-Apcr+] 

e_T0( } = < e ~ T > = jL ^T3 FT1 f 7H\ ( 26 ) 

so that, recalling that ((1) = by construction, 

0(A) =-max[C(p)-Apa+] . (27) 



The inversion of the Legendre transform yields 



C(p) - nun [A P( t + - 0(A)] , (28) 



and it is easy to check that the fluctuation relation is equivalent to 0(A) = 0(1 — A). 



C. Internal symmetries and the fluctuation relation 

Assume that there exists a map I on the space of trajectories r(t) such that I 2 — 1 and that the measure Vr is 
inviariant under /, i.e. Vlr = Vr. Then, consider a segment of trajectory r(t), t e [— r/2,r/2] and define 

S in^M (29) 

where V[r(t)} is the probability of observing r(t), in the stationary state, for t € [— r/2,r/2] irrespectively of what 
happens outside the interval [— r/2, r/2]. It is easy to show that the PDF of S T verifies the fluctuation theorem. Indeed 

(e~ XS -) = J VrV[r(t)}e- xs - = J VrV[r(t)] 1 - x V[Ir(t)] x = J VrV[Ir(t)] 1 - x V[r(t)] x = (er^-^) . (30) 

Thus, if the limit 

0(A) = - lim t- 1 \n(e- xs ^) (31) 

exists, it verifies the relation 0(A) = 0(1 — A) from which the fluctuation relation for the PDF of S T follows. Lebowitz 
and Spohn showed that the limit 0(A) indeed exists for generic Markov processes and it is a concave function of A. 
Moreover they showed that S T can be identified with the entropy production rate -over the time interval r- in the 
stationary state up to boundary terms, i.e. terms that do not grow with t, if / is chosen to be the time reversal, 
Ir(t) = r(-t). 



D. Entropy production rate 



We are interested in the explicit form of S T for the equation of motion (11) in the case in which g a fi{t) = S a pg(t) 
and h(t) = h[f(t)] is an external nonconservative force that does not explicitly depend on time: e.g., in d = 2, 
h = a(—y,x). Note that the functions v(t) and g(t) are such that v(t) = v(—t) while g(t) is proportional to 0(t), and 
both decay sufficiently rapidly in time. The probability distribution of the noise p(t) is 



V[p(t)] oc exp 



-i J dtdt' p a {t)v~\t-l!)p a {l!) 



(32) 



where v 1 (t) is the operator inverse of v{t), see Eq. (20). The probability distribution of f(t) is obtained substituting 
p(t) obtained from Eq. (11) in Eq. (32). One has 



V[f(t)} cxexp j-i J dtdt' mr a (t) + J dt" g{t - t")r a {t") + ^W)] - K[r{t)\ 

■ a (t>) + J dt"'g(t' - t"')r a (t"') + ^rW)] - KW)] 



(33) 
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After some algebra it is easy to see that 



V[r{— t)] ocexp 



dtdt' 



mr a (t) - J dt"g(t" - t)r a (t") + - KW)] 

a (t') - J dt"'g(t"' - t')f a {t"') + ^rW)\ - h *W)\ 



(34) 



To compute S T we should consider the probability of a segment of trajectory [— r/2, r/2] and then send r to oo, 
neglecting all boundary terms. As the functions g(t) and v(t) have short range, the trajectories f(t) decorrelate, 
say, exponentially fast in time and up to boundary contributions one can simply truncate the integrals in V[r(t)] in 
t,t'e [-r/2, r/2]. 

Let us now discuss the contributions to — \nV[f(— t)] + In V[f(t)] that do not trivially vanish. One has: 



a "kinetic" term of the form 



J dtdt' mf a (t)i/ _1 (t-t / ) J dt" 'g(t" - t')r a {t") + mr a (t)i/ _1 (t - t') J dt" g(t' - t")r a (t") 
= J dtdt' mr wit)^ 1 (t-t') J dt" f(t' -t")f a {t") = J dtdt'mf a {t)T- 1 (t-t')r a {t') . 



(35) 



If the bath is in equilibrium, this term trivially vanishes as it is the integral of the total derivative of the kinetic 
energy. But it also vanishes for a nonequilibrium bath. Indeed, by integrating by parts first in t and then in t', 
we find 



J dtdt'r a (t)T- 1 (t-t')r a (t') = - J dtdt' r^jT^it - t')r a (t') = 

J dtdt' r a {t)^T-\t - t')r a {t') = - J dtdt' r a {t)T-\t - t')r a (t') = . 

where we used that T _1 (t) is even and short ranged and we neglected boundary terms, 
a "friction" term of the form 

\ J dtdt'dt"dt"' [f a (t")g(t" - t)v~\t - t')g(t"' - t')r a (t"') - r a (t")g(t - t") v -\t - t')g(t' - t"')f a (t"')} 
This term vanishes because the function 



(36) 



K(t" - t'") = J dtdt' g{t" - t)v-\t - t')g(t"' - t') 



(37) 
(38) 



is even in its argument as one can easily check. 
• a "potential" term of the form 

pt/2 f-r/2 f-r/2 



rT/Z rT/Z pT/'Z fly 

= - / dt dt' dt" f{t-t")r a {t")v-\t-t')^-[r{t')\ 

J— r/2 J-t/2 J-t/2 ar a 

/ T / 2 r T l 2 8v 
dt dt'T-\t-t')f a (t)^-[f(t')] . 
-r/2 J-t/2 O r a 



(39) 



This term is related to the work of the conservative forces. If the bath is in equilibrium, it vanishes being 
the total derivative of the potential energy. It vanishes also for a harmonic potential V(r) — \kr 2 because 
■§^-[r(t)] = kr a (t) and one can use the same trick used in Eq. (36). It docs not vanish in general. 

a "dissipative" term 



f-r/2 r-T/2 

Sf ss = dt dt'T-^t-t'y^Klrit')}. 

J-t/2 J-t/2 



r/2 
-r/2 



(40) 



This term is related to the work of the dissipative forces. If the bath is in equilibrium at temperature T, this is 
exactly the work of the dissipative forces divided by the temperature of the bath. Otherwise, the work done at 
frequency ui is weighted by the effective temperature at the same frequency. 
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The expression of the total entropy production over the interval [— t/2, t/2] is then 



f-r/2 r-T/2 

S T = S^+Sf ss = / dt / dt'T-\t-t')r a {t) 

J-t/2 J-t/2 



/2 
/2 



dV , 



dr n 



f(t')] + h a [r(t')] 



f-r/2 f -r/2 

= dt dt'T- l {t-t')r a {t)F a {t') , 

J-t/2 J-t/2 



-t/2 J-t/2 



(41) 



where F a (t) = h a [r(t)] — ^-[f{t)} is the total deterministic force acting on the particle at time t. 
The latter expression can be rewritten as 



f-r/2 r-r/2 

S T ~ dta t = dt (aY + af ss 

J-t/2 J-t/2 



/2 
t/2 



(42) 



defining a entropy production rates a t , aY and af lss (modulo subdominant terms in the large t limit) as 



~- ~V i „diss 
°t = O t + t 



f dt'T-\t-t')[f a {t)F a {t')+r a (t')F a {t)] , 
^ = dt' T~\t - t') [f a (t)-^-[r(t')] + r a (t')^-{r(t] 

af ss = f di'T- 1 (i-i')[ra(t)/ia[r(i')]+' i a(0^[^)]] • 

J — OO 



(43) 



We recall that if the bath is in equilibrium this expression reduces to the work done by the nonconservative forces 
divided by the temperature of the bath, as expected. If the bath is not in equilibrium, but the potential is harmonic, 
only the contribution a^ lss of the nonconservative force has to be taken into account. The reason why the work 
of the conservative forces produces entropy if the bath is out of equilibrium and the interaction is nonlinear is that 
the nonlinear interaction couples modes of different frequency which are at different temperature, thus producing an 
energy flow between these modes; this energy flow is related to the entropy production. 

It is also important to remark that boundary terms, that are usually neglected, can have dramatic effects on the 
large fluctuations of S T even for r — > oo, as pointed out by Van Zon and Cohen [42]. This happens if the PDF of 
the boundary term has exponential tails. Thus, boundary terms cannot be always neglected, at least for very large 
values of p. A good empirical prescription to remove boundary contributions is the following: in equilibrium a t must 
be a total derivative as no dissipation is present. So, removing a total derivative (a boundary term) one can define 
S T in such a way that it vanishes identically in equilibrium. This definition turns out to be the one that verifies the 
fluctuation relation for all \p\ < p* , p* being the maximum allowed value of p for t — > oo [42, 43]. We shall discuss 
this point in more detail later. 



III. A DRIVEN PARTICLE IN A HARMONIC POTENTIAL 



In this and the next Section we discuss some examples on which we test the fluctuation relation for S T defined 
in (41). We first consider the simplest fully analytically solvable case in which there are no applied forces and the 
potential is quadratic. We derive the fluctuation-dissipation relation between induced and spontaneous fluctuations 
in the position of the particle and we relate it to the time-dependent temperature of the bath defined in Eq. (19). We 
then show that in this simple case the fluctuation relation for S T in (41) reduces to the usual one with the temperature 
of the bath. 



A. The fluctuation-dissipation relation 



In the harmonic Brownian particle problem with no other applied external forces the dynamics of different spatial 
components are not coupled. Thus, without loss of generality, we henceforth focus on d = 1. In Fourier space, the 
Langevin equation reads 



-mti) 2 i(aj) — iojg(Lu)x(ui) = —kx(ui) + p(uj) 



(44) 
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with the noise-noise correlation 

( p(lu) P (lu') ) = 2ttS(uj + u')u((j) . (45) 

The linear equation (44) is solved by 

x(u) = G(w)p(w) , G(w) = 2 1 , r— r , (46) 

—into 1 — iLug(LU) + k 

and one finds the correlations 

{x{u))x{u')) = G{w)G{-uj)2-k5{u) + J)v{u>) , 

( x(u)p{u') ) = G(lo)2t:8(lu + w>(w) . (4?) 
Note that G(w)G(-w) = |G(u;)| 2 ; then 

( .x(w)a;(w') ) = C(w)2tt(5(w + w') with C{w) = \G{uj)\ 2 v{uj) . (48) 

In a problem solved by 

/oo 
dt' G(t - t') [p(t') + h(t')} + IC , (49) 
-oo 

adent linear response is 

G(t - t') , (50) 



where IC are terms related to the initial conditions, the time-dependent linear response is 

6(x(t)) 



R(t - 1 1 ) 
and 



Sh(f) h=0 



/OO 
dt e lut R{t) = G(w) 
-oo 



R(u) = / dt e tuJt R(t) = G(uj) . (51) 



Note that the response function is related to the correlation (x(t)p(t') ) by Eq. (47): 

2ttS(lu + u')R{uj)v{uj) = ( x(lu) P (lu') ) . (52) 

Now, we can check under which conditions on the characteristics of the bath [g(t — t') and v(t — t')] the fluctuation- 
dissipation theorem (for the particle) holds and, when it does not hold, which is the generic form that the relation 
between the linear response and correlation might take in this simple quadratic model. Eq. (46) implies [40] 

hnR(w) = ImG(w) = to Reg(uj) \G{uo)\ 2 , (53) 

and then using equation (48) 

loC{lo) v(lo) 



2lmR(uj) 2Rc.gH ^ ' ^ 54 ^ 

We see that the fluctuation-dissipation theorem holds only if this ratio is equal to T, see Eq. (16). Otherwise, the 
relation between linear response and correlation of the particle is given by the frequency-dependent temperature of 
the bath, T(oj), defined in (18). The measure of the modification of the fluctuation-dissipation theorem given in (54) 
is the effective temperature of the system. The use of this name has been justified within a number of models with slow 
dynamics and a separation of time-scales [9, 10] but it might not hold in complete generality [11]. It is important to 
remark that the effective temperature in the frequency domain, T(u>) is not equal in general to the Fourier transform 
of the effective temperature T(t) -defined in the introduction- which is the ratio between the noise-noise correlation 
and memory function in the time domain. 

Let us now discuss some environments that we shall use in the rest of the paper. 



1. Equilibrated environments 

For any environment such that the right-hand-side in Eq. (54) equals T the fluctuation-dissipation theorem holds. 
In the time domain, this condition reads 

Tg(t) - 9(t)v{t) , v{t) = T[g{t) + g(-t)} = Tg(\t\) . (55) 

In particular, this is satisfied by a white noise for which v{t) = 2T-fS(t) and g(t) = 2-fS(t)9(t) (remember that 
0(0) = 1/2). The fluctuation-dissipation theorem also holds for any colored noise - with a retarded memory kernel g 
and noise-noise correlation v - such that the ratio between Keg(oj) and v(oj) equals (2T) _1 . This requirement applies 
to any equilibrated bath. 
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2. Nonequilibrium environments 

Instead, for any other generic environment, the left-hand-side in Eq. (54) yields a non-trivial and, in general model- 
dependent, result for the effective temperature. 

A special case that we shall study in Appendix B is the one of an ensemble of N equilibrated baths with different 
characteristic times and at different temperatures. In this case, the noise p in Eq. (11) is the sum of N independent 
noises, 

N 

P = ^Pi , (Pia(t)Pjp(t')) = 5 a p6ijTiVi{t-t') , (56) 

i=l 

and the friction kernel is given by 

s(t-f) = I>(t-f) ■ ( 5? ) 



We have extracted the temperature Tj from the definition of Vi(t) in order to simplify several expressions. As the pi 
are independent Gaussian variables, p = J2i Pi i s still a Gaussian variable with zero mean and correlation 

(p a (t) P 0(t')) = 5 af3 ]T T lVl (t - t') . (58) 

i 

Thus, in the Gaussian case the N equilibrated baths are equivalent to a single nonequilibrium bath with correlation 
given by Eq. (58) and friction kernel given by Eq. (57). In frequency space we have 

N N 

ffM = ^2di(u) , v{uj) = ^TiVi{u) , (59) 

i=l i=l 

with 

Vi(w) = 2Rc 5l (w) , (60) 
as each bath is equilibrated at temperature Tj. The frequency-dependent temperature is then given by 

TM = • (61) 

Note that if the functions Vi(u)) are chosen to be peaked around a frequency Ui, by choosing suitable values of and 
T, one can approximate a single nonequilibrium bath with N baths equilibrated at different temperatures. 



B. Large deviation function 

We now compute the large deviation function in the harmonic case, V(r) = ^kr 2 . In this case uX is a total 
derivative and only the term cr* ss , related to the nonconservative forces, is relevant. We show that the characteristic 
function 0(A) of af lss exists, is a convex function of A and verifies the fluctuation relation 0(A) = 0(1 — A). 



1. Equilibrium bath 

As a first illustrative example we consider the case of an equilibrium white bath. The model we study is a 
two dimensional harmonic oscillator with potential energy V(x,y) — |(x 2 + y 2 ) coupled to a simple white bath in 

equilibrium at temperature T, and driven out of equilibrium by the nonconservative force h = a(—y, x). The equations 
of motion are 



m'x t + jx t = -kx t - ay t + & , 
myt + jy t = -ky t + ax t + n t , 



(62) 
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where £t and 77* are independent Gaussian white noises with variance (£t£o) = (?7t?7o) = 2jTS(t). The memory 
friction kernels g a p(t — s) are simply 5 a pg(t — s) = 25 a pjS(t — s)6(t — s) in this case, with 7 the friction coefficient. 

Defining the complex variable z t — (x t + iyt)j\f2 and the noise p t — (£t + iVt)/V2 the equations of motion can be 
written as 



mz t + "fz t = —KZ t + pt , (63) 
where n = k — ia, (ptPo) = (PtPo) = an d (PtPo) = 2-fTS(t). The complex noise p t has a Gaussian PDF: 



V[pt] oc exp 



1 r 



dt p t pt 



The energy of the oscillator is H = mz'z + kzz, and its time derivative is given by 

dH 

-^-j- = 2mRez t z t + 2kRez t z t = 2almz t z t - 2jz t z~ t + 2Rei t p t =W t -W t 



(64) 



(65) 



where Wt — 2almi t z t = a(x t yt — yt%t) is the power injected by the driving force and Wt = 2jitz t — 2Rez t p~t is the 
power extracted by the thermostat (henceforth we choose the sign of the power in order to have positive average). 

The entropy production rate (43) reduces, as expected, to the injected power divided by the temperature, at = (3Wt, 
where [3 — 1/T (one could also consider the entropy production of the bath, a t — /3Wj; for completeness we discuss 
it in Appendix C). 

We want to compute the probability distribution function (pdf) of the entropy production rate u t = (3W t - The 
average value of at is in this case given by a + = 2a 2 / (-fk). From Eq. (41) we can rewrite the total entropy production 
S T in terms of the complex variable z t : 



rr/2 

/ dt a t = 2a(3 Im 



-r/2 



r/2 



dt z t z t ■ 



r/2 



(66) 



As already discussed, it is easier to compute the characteristic function 4>(X), Eq. (31), in terms of which the fluctuation 
relation reads </>(A) = <p(l — A). To leading order in r we can neglect all the boundary terms in the integrals. After 
integrating by parts we have 



72 



S T = 2af3i I dt z t z t ■ 

J -t/2 

In terms of the PDF of the noise (64) we obtain 

(exp[-A5 r ]) =M~ 1 j d V[pt] exp 



2ia\ [ T/2 

T/2 



dt z t zt 



(67) 



(68) 



and the normalization factor J\f = J dV[p t ] is simply given by the numerator calculated at A = 0. 

To leading order in r the function <^>(A) should not depend on the boundary conditions in Eq. (68). Thus, we impose 
periodic boundary conditions, z(t/2) = z(—t/2) and i(r/2) = z(— t/2), and we expand Zt in a Fourier series, 



where Auj — 2tt/t and u n = nAui. For r — > 00 

dio 
,2^ 



2tt ^ 



Z 7r 6 



-lUJnt 



Zt 



e~ luJt z 



dt e l0Jt z t 



and the equations of motion become 



Pu 



_ Pu 



-uj 2 m + k — iuj"f D(uj) 



(69) 



(70) 



(71) 
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Note that in the limit a = the Green function G(a,uj) — 1/D(u>) reduces to the one used above to compute the 
violation of the fluctuation-dissipation theorem induced by a nonequilibrium bath. The distribution of the noise is 
given by 



V[p ul } = cxp 



2 7 r J-, 



du> 
2^ 



exp 



1 Alu ^ 
2^T 27 ^ PnPn 



(72) 



Substituting Eqs. (69) and (72) into Eq. (68) we get 



(exp[— XS T ]) = N 1 J dp n exp 
and using Eq. (25) 



Alu 
"~2tT 



2-^ \ 9-t/T 



2aXm n \p n \ 2 
V2 7 T T\D{uu n )\ 2 



n 



\D(u; n ) 



(73) 



To show that (j) verifies <p(X) 



lim 

T— >CC 



4 7 aAw„ 



|£>(w„ 



do; 

— In 

2tt 



47a Aw 



|D( W ) 



(1 — A) and hence the fluctuation theorem, note that 



„(A)-r,(l-A)= / ^ln 



\D{uu)\ 2 - A^aXuu 
\D{lo)\ 2 - 4 7 a(l - X)lu 



(74) 



(75) 



and, as |Z?(w)| 2 — 40:70; = \D{— lu)\ 2 , the integrand is an odd function of to and the integral vanishes by symmetry. In 
Appendix A we show that the same result is obtained if one uses Dirichlet boundary conditions (at least for m = 0, 
where the computation is feasible); this result supports the approximations made when neglecting all the boundary 
terms in the exponential. Moreover, in the case m = the large deviation function ((p) can be explicitly calculated; 
defining t — -f/k and er = a + r /2 = a 2 /k 2 , we obtain 



C(P) = r 

Thus, for t — > 00 the PDF of p has the form 



1+P<t - V(l + ct )(1 +P 2 c ) 



7T T (p) (X CXp 



T - 

— C(P)^o) 
TO 



(76) 



(77) 



Note that To is the decay time of the correlation function of z t [i.e. (ztzo) oc exp(— t/ro)] and do is the average entropy 
production over a time ro/2. Thus, r is the natural time unit of the problem (as expected); remarkably, the function 
C = t C depends only on a and not on the details of the model. It would be interesting to see whether the same 
scaling holds for more realistic models. 

In summary, we found that for all driving forces, i.e. all values of a, the fluctuation theorem holds for the entropy 
production rate (41). For a white equilibrium bath this result has already been obtained in [6]. The temperature 
entering the fluctuation theorem is the one of the equilibrated environment with which the system is in contact, 
although it is not in equilibrium with it, when the force is applied. 

Let us also stress that one can easily check that the fluctuation-dissipation relation holds in the absence of the drive 
(see Sect. Ill A) but it is strongly violated when the system is taken out of equilibrium by the external force. 



2. N on- equilibrium bath 



We now generalize the calculation to the case of a generic nonequilibrium bath; the equation of motion becomes 

/oc 
dt' g(t - t')z v = -KZ t + p t , (78) 
-00 

where as before n — k — ia and (ptPo) = v (t)- The functions u(t) and g(t) are now arbitrary (apart from the condition 
g(t) = for t < 0), hence they do not satisfy, in general, Eq. (55). Note that Eq. (78) provides a model for the 
dynamics of a confined Brownian particle in an out of equilibrium medium [35-37]. 
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The dissipated power is given by 

dH n T . _ 
—— = 2a lmz t z t 
at 



2Rc 



/OO 
dt' 



$(i - t')z t z v + 2Rez t p t =W t -Wi 



t > 



(79) 



where as in the previous case W t — 2almz t z t is the power injected by the external force and W t = 2Re dt' g(t 

t')z t Zt> — 2Rcz t pt is the power extracted by the bath. 

For the harmonic model is a boundary term and Eq. (41) gives 



9tt £—1 



27T ^ T{u) n ) 

Tl — — OO V ' 



,r/2 r-T/2 

2ai / dt dt' T~ x {t - t')z t z t > . 

J-t/2 J-t/2 



(80) 



Note that the last equality holds neglecting boundary terms. 

Let us now compute <j>dissW- The computation is straightforward following the strategy of section III. In Fourier 
space, Eq. (78) reads 



Pu 



The probability distribution of p u is 



Thus [see Eqs. (72) and (73)], 



(exp[-XSf ss }) ^Af- 1 I d Pn exp 



muj 2 + k — iujg(uj) D(u>) 
duj \ PuJ \ 2 ' 



V[pJ\= cxp 



/ 

J —i 



2n v(lo) 



(81) 



(82) 



Alo 



OO 

E 

i— — OO 



\Pn\ 



2a\uj n \p n \ 2 
v(u n ) T(u n )\D(uj n )\ 2 



n 



2a\uj n v(uu n ) 
T(uj n )\D(u n )\' 



and using the definition of T{u) given by Eq. (18) 



(A)= / 



duo 
2n 



In 



1 - 



4a\ujReg(uj) 



\D(uW 



(83) 



(84) 



It is easy to prove that |£>(a->)| 2 — 4au>Reg(uj) = \D{— oj)\ 2 . Using now the same trick employed in Eq. (75), one shows 
that fidissW = 0dj SS (l - A). 

An alternative definition of entropy production rate in which one assumes that it is proportional to the power 
injected by the external drive, af = Q~ 1 W t , via a parameter 6 which has the dimension of a temperature, has been 
often used in the literature [22-24, 44]. With this definition, the total entropy production over a time r is given by 
(neglecting boundary terms) 



2ai 



72 J • „ Au; 
at z t z t — —2a- 

r/2 



2tt 



E 



e 



(85) 



i.e. T(uj) is replaced by a constant 9 that is taken as a free parameter that one adjusts in such a way that the PDF of 
Sf is as a close as possible to verify a fluctuation relation [22-24]. Substituting T(uj) with a constant <d in Eq. (83) 
one obtains 



2a\ujv(u)) 

e|£H| 2 



(86) 



However, it is not possible to find a value of such that 4>q(X) satisfies the fluctuation theorem for the harmonic 
problem in contact with a generic bath. We shall show in section IV B that the use of a single parameter 6 constitutes 
a rather good approximation when the dynamics of the particle occurs on a single time scale. 

In conclusion, the fluctuation theorem is satisfied when the entropy production rate is defined using the power 
injected by the external drive with the temperature of the environment defined as in (18). 

In Sect. Ill A 2 we also introduced a complex bath made of many equilibrated baths at different temperatures, 
eventually acting on different time scales. In Appendix B we prove that, as expected, the PDF of 5* ss defined in 
Eq. (80) also verifies the fluctuation theorem in this case - while the PDF of S® does not. For such a multiple bath 
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one can also consider the entropy production of the baths, defined as the power extracted by each bath divided by 
the corresponding temperature. This quantity is of interest if one could identify the different thermal baths with 
which the system is in contact; clearly, this is not possible in glassy systems where the effective temperature is self- 
generated. Nevertheless, the study of systems of particles coupled to many baths at different temperature is of interest 
in the study of heat conduction. In Appendix C we prove that the entropy production rate of the baths verifies the 
fluctuation theorem, at least for \p\ < 1 (see also [42, 43]). 



IV. A DRIVEN PARTICLE IN AN ANHARMONIC POTENTIAL: NUMERICAL RESULTS 

In this section we investigate numerically Eq. (78) for a particular choice of the nonequilibrium bath and in presence 
of a linear and nonlinear interaction. In the linear case, we find that the numerical results reproduce the analytical 
results of the previous section. This finding confirms that the boundary terms we neglected in the analytical compu- 
tation are indeed irrelevant. In the nonlinear case, we find that the fluctuation relation seems to be satisfied quite 
well for S^ lss , although strictly speaking it only holds for S^ tss + . We shall discuss the reason for this below. 

We consider the simplest non trivial case, in which a massless Brownian particle is coupled to two equilibrated baths: 
a white (or fast) bath at temperature Tf and a colored (or slow) bath with exponential correlat ion at temperature T s . 
This model has been studied in detail in [10] and is relevant for the description of glassy dynamics when the time scales 
of the two baths are well separated, as will be discussed in Sect. V. The equations of motion are given by Eq. (78) with 
9(t) = 9f(t) +9 s (t), 9f{t) = -f f S(t) and g s (t) = 0(t)^e~^ , or, equivalently, g f (u) = j f and g s (u) = j s /(l - iivr s ). 

We use a generic rotationally invariant potential V(x,y) = V % v = V(\ z \ 2 )- The noise is the sum of a fast and 
a slow component. Then Eq. (78) becomes 

/t , 
dt' e-^Tzv = -z t V{\zt\ 2 ) + iaz t + p{ + p\ , (87) 
-oo 

where (p{ p[,) = 2-ffTfS(t - t'), {p s t p s t ,) — 2kr» e -l*-*'l/^ anc i V'(x) is the derivative of V(x) wtih respect to x. It is 
convenient to rewrite this equation as 



-if z t = -z t V'{\z t \ 2 ) + iaz t + p[ + u t 



_ IsZt 



where we introduced the auxiliary variable u t and v t is a white noise with correlation (v t v t >) — 2j s T s 5(t — t'). The 
power injected by the external force is, as usual, Wt = 2almi t z t , while the power extracted by the two baths can be 



written as w/ = 2Re 



zt yis'zt- p{ 



and W t s = 2Re 



z t -trz t - u t 



For concreteness we focus on the potential V(\z\ ) = f|z| and compare with the results obtained for the harmonic 
case, V(|z| 2 ) = fc|z| 2 . The simulation has been performed for a = 0.5, Tf = 0.6, 7/ = 1, T s = 2, 7 S = 1 and t s = 1; 
we set k = 1 in the linear case and g = 1 in the nonlinear one. The system (88) is numerically solved via a standard 
discretization of the equations with time step St = 0.01; the noises are extracted using the routine GASDEV of the C 
numerical recipes [45]. 

We found that and S^ lss are uncorrelated (within the precision of the numerical data). To the extent that this is 
the case, their PDFs can be studied separately. Unfortunately, the PDF of is too noisy to allow for a verification of 
the fluctuation relation in the nonlinear case. This is probably due to the fact that while in the linear case reduces 
to a boundary term [46] , in the non- linear case 'spurious' boundary contributions might be difficult to eliminate [42] . 
Indeed, for the accessible values of r, the variance of is much larger than its average (while the fluctuation relation 
predicts a variance of the order of a\). The lar ge variance might be a finite-r effect due to a boundary term with 
fluctuations contributing to the ones of but not to the average. If this were the case, the fluctation relation should 
hold for \p\ < 1 and very large r. However, the required values of t might be so large to render the fluctuation relation 
unobservable in practice, see [20, 43]. For this reason, we shall not discuss the data for . The validity of the FR 
for (possibly minus a boundary term) in the nonlinear case remains an open question that should be addressed 
by future work. 
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FIG. 1: PDF of of ss : a) The large deviation function for the harmonic potential at r = 10, 20, 30; the full line is the analytical 
result, b) The function f(p) = [C,diss(p) — Cdiss(— p)]/o+ lss for the harmonic potential: the fluctuation theorem predicts a 
straight line with slope 1, represented by a full line, c) The large deviation function for the quartic potential at r = 10, 20, 30. 
d) The function f(p) for the quartic potential: also in this case the fluctuation theorem is well verified. 



A. Entropy production rate 

Let us now discuss the behavior of Sf ss . The dissipative contribution to the entropy production rate, see Eq. (41), 
is given by 



diss 
~t 



I dt 1 T^(t - t') 2aha[z t z t ' + z t >z t ] • (89) 

J — OC 



The inverse of the frequency dependent temperature, l/T(o;), is 

1 lf {l+UJ 2 T 2 s )+ ls 



T(u>) T nf (l + LJ*T*) + T sls ' 

see Eq. (61). Thus 



(90) 



T-Ht) = ±6{t)+*(l-*!) e ^-, with n=' Tnf+TslS > (91) 
T f T HS T s \ T f J 20 t s y 7/7/ 

and T _1 (i) decays exponentially for large t. Note that, if the bath is in equilibrium, T s — Tf = T, one has 
T _1 (t) = 5{t)/T and of ss = 2aknz t z t /T = W t /T as expected [recall that in our convention J**^ dt' 5(t - t') = \]. 

The data for Sf ss are shown in Fig. 1. The large deviation function Cdiss(p) is reported in panel a) for the harmonic 
and in panel c) for the quartic potential. The average af ss is equal to 0.332 in the harmonic case and to 0.276 in the 
quartic case. The function Cdiss(p) converges fast to its asymptotic limit r — * oo (note that even the data for r ~ 10 
are in quite good agreement with the analytic prediction for the harmonic case). The fluctuation theorem predicts 
f(p) = {Qiss(p) — (diss(—p)]/of ss = p. The function f(p) is reported in panel b) for the harmonic and in panel d) 
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FIG. 2: PDF of of: a) The large deviation function for the harmonic potential at r = 10, 20,30; the full line is the analytical 

result, b) The function f(p) = [C(p) — C^f 1 )]/ "^ f° r the harmonic potential. The full line is the analytical prediction, the 
dashed line is the fluctuation relation, the dot-dashed line has slope Tf/T s . c) The large deviation function for the quartic 
potential at r = 10, 20, 30. d) The function f(p) for the quartic potential; the dashed line is the fluctuation theorem, the 
dot-dashed line has slope Tf/T s . 



for the quartic potential. In the harmonic case the numerical data are compatible with the validity of the fluctuation 
theorem, as predicted analytically. Remarkably, the same happens in the quartic case for which we do not have an 
analytical prediction. 

These results support the conjecture that, if S^ ss and are uncorrelated, the PDF of S^ zss verifies the fluctuation 
theorem independently of the form of the potential V(x, y). 



B. Approximate entropy production rate 

We also investigated numerically the fluctuations of the entropy production Sf used in some numerical and exper- 
imental studies [22-24, 44] and that we discussed in section IIIB 2. For this model it is given by 

*? = f = |W ■ (92) 

Rather arbitrarily we set = Tf in the definition of erf . This reflects what is usually done in numerical simulations, 
where the dissipated power is divided by the "kinetic" temperature, i.e. the temperature of the fast degrees of freedom. 
Note that the choice = Tf does not affect the function Ce(p) since the variable p is normalized [i.e., Ce(p) = C(p) 
does not depend on 0, see Eq. (1)] but it changes the average that is proportional to _1 . 

The data for af are reported in Fig. 2. The harmonic case is shown in panels a) and b) while the anharmonic case 
is presented in panels c) and d). We have ctv' = 0.455 for the harmonic potential and av f = 0.366 for the quartic 
one. The numerical result for the large deviation function of erf agrees very well with the analytical prediction in the 
harmonic case but, as discussed in section IIIB 2, it does not verify the fluctuation theorem for = Tf, as one can 
clearly see from the right panels in Fig. 2. 
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t 

FIG. 3: Normalized autocorrelation functions of zt for the harmonic oscillator with a = 0.5, k = 1, Tf = 0.6, T s = 2, 7/ = 1, 
7 S = kr a and t s — 1, 50, 250. 

Remarkably, in both the harmonic and anharmonic cases the function f(p) = [C{p) ~ C(~p)]/ (T + i is approximately 
linear in p with a slope X such that 1 > X > Tf/T s , i.e. — £(— p) ~ X pcrT^ . If /(p) ~ Xp, one can tune the 
value of in order to obtain the fluctuation relation — £(— p) = per®, by simply choosing = Q e ff = Tf /X, thus 
defining a single "effective temperature" e // <E [T/,T S ]. From the data reported in Fig. 2 we get a slope X ~ 0.66, 
that gives 9 e// = T//X - 0.9. 

This behavior reflects the one found in some recent experiments [22-24, 41] in situations in which the dynamics of 
the system happens essentially on a single time scale. This is the case also in our numerical simulation: in Fig. 3 we 
report the autocorrelation function C(t) = Re(ztZo) of z t (computed in Appendix D) for the harmonic potential. The 
present simulation refers to the curve with r s = 1, which clearly decays on a single time scale. 

In Fig. 4 we report the parametric plot x(C) ( see t ne Introduction and Sect. Ill) for the same set of parameters, 
but a = 0. The integrated response is given by = J* dt' R(t') and R(t) is computed in Appendix D. We see that, 
for t s = 1, the function x(C) nas a slope close to —1/Tf at short times (corresponding to \ ~ 0). For longer times, 
the slope moves continuously toward —1/T e j /, with T e f f ~ 1.37. This value of T e f f is of the order of 7/ ^^ sTs =1.3, 
which means that on time scales of the order of the (unique) relaxation time the two baths behave like a single bath 
equilibrated at an intermediate temperature. This would be exact if the time scales of the two baths were exactly 
equal. 

It is worth to note that in this situation we get T e ft ^ ©e//, that is, the effective temperature that one would 
extract from the approximate fluctuation relation of Fig. 2 does not coincide with the effective temperature obtained 
from the x(C) plot of Fig. 4. In particular, we get Tf < © e // < T e ff. this relation is consistent with the results of 
[41] obtained from the numerical simulation of a sheared Lennard- Jones-like mixture, even if the coincidence might 
be accidental. 



C. Discussion 

Let us summarize the results in this section. The numerical simulation of the non-linear problem confirms that the 
fluctuation theorem is satisfied exactly when the entropy production rate af lss is defined using the power injected by 
the external drive and the temperature in (18) is used. 

In situations in which the dynamics of the system happens on a single time scale, a fitting parameter O can be 
introduced to obtain an approximate fluctuation relation on the entropy production rate Wt/Q. However, is not 
necessarily related to the effective temperature T e f f that enters the modified fluctuation-dissipation relation. In the 
systems considered so far [41] one finds < T e ff. However, this is just an approximation that fails in more generic 
non-equilibrium situations. In the next section we show that, when the dynamics happens on different, well separated, 
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FIG. 4: Parametric plot of the integrated response \if) as a function of the correlation function C(t) for the harmonic oscillator 
with a = 0, k = 1, Tf = 0.6, T s = 2, 7/ = 1, 7 S = fcr s and r s = 1 (continuous curve) and t 3 — 250 (broken curve). The 
dot-dashed line has slope —1/1.37, the dashed lines have slope — 1/T S and — 1/T/. 

time scales, it is impossible to find a single value of O such that af — W t /Q verifies the fluctuation relation. 

V. SEPARATION OF TIME SCALES AND DRIVEN GLASSY SYSTEMS 

In this Section we discuss the application of our results to glassy systems. After presenting the general argument, 
we show explicitly that the dissipative entropy production satisfies the fluctuation relation for the p-spin spherical 
model. We finally discuss an adiabatic approximation that allows one to derive approximate results in the case of 
systems with well-separated time-scales. 

A. Background 

As discussed in the Introduction, in the study of mean-field models for glassy dynamics [10, 34] and when using 
resummation techniques within a perturbative approach to microscopic glassy models with no disorder, effective 
equations of motion of the form of Eq. (10) are obtained: 

jcpt = -fi(t) (fn+ f dt 1 E(t, t') tp v +p t + ah[(p t ] , (93) 

where p t is a Gaussian noise such that (ptpt'} = D(t,t'). The "self-energy" E(£, t') and the "vertex" D(t,t') depend 
on the interactions in the system throught the correlation and response of the field (p. In absence of drive (a = 0) 
this equation has a dynamic transition at Td separating a high temperature phase in which the dynamics rapidly 
equilibrates from a low temperature phase in which the dynamics is non-stationary and the fluctuation-dissipation 
relation is violated. This means that the time needed for the system to equilibrate is the longest time scale that is 
unreachable in the calculation (it already diverged with TV) . 

In the case of a driven mean field system [26-28] , the external force is also present in Eq. (93) and after a transient 
the system becomes stationary for any temperature, i.e. p{t) = p, E(i, t') = E(i — t'), and D(t, t') = D(t — t'). The 
functions D and E depend on the strength a of the driving force. In this case, Eq. (93) resembles Eqs. (78) and (87), 
and our results of Sect. Ill apply. By analogy with Eq. (18), the effective temperature is defined in terms of E and D, 
see below. 

As discussed in [26-30], for small a the system shows a completely different behavior above and below Td, reflecting 
the presence of a dynamical transition at a = 0. Above Td, the fluctuation-dissipation theorem holds in the limit 
of a — > 0; the transport coefficient related to the driving force a approaches a constant value for a — > (the linear 
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response holds close to equilibrium) and the systems behaves like a "Newtonian fluid". Below Td, the fluctuation- 
dissipation relation is violated also in the limit a — -> and the transport coefficient diverges in this limit: the system 
is strongly nonlinear. For a wide class of systems, see the Introduction and Ref. [34], the relation between D and £ 
takes a very simple form in the limit a — > 0: the effective temperature T e ff(t — t') defined from the ratio between 
induced integrated response and correlation (see the Introduction) is given by the temperature of the bath T for small 
1 1 — t' | and by a constant T e // > T for large 1 t — t' |. Thus, we expect that above Td (and for a <~ 0) the system 
behaves as if coupled to a single equilibrium bath (and the fluctuation theorem holds for the entropy production rate 
defined as a t — ah[<pt]<pt/T), while below Td the system behaves as if coupled to two baths acting on different time 
scales and equilibrated at different temperatures. 

As already remarked in the introduction, these single-spin equations of motion are valid for times r that are finite 
with respect to the size N. For these times the system behaves like independent spins moving in a harmonic potential 
in contact with a nonequilibrium environment. Thus, for S T with r in this regime the results of section III apply and 
the correct definition of the entropy production rate is given by Eq. (80), i.e. by the power injected by the external 
force alone, divided by the frequency-dependent effective temperature. 



B. The spherical p-spin 

The (modified) fluctuation relation can be checked explicitly for the p-spin spherical model. This model realizes 
explicitly the behaviour described above. The asymptotic dynamics in the low temperature phase occurs in a region 
of phase space that is called the threshold and it is far from the equilibrium states [3] since times that grow with TV 
are needed to reach them. 

The effective equations of motion for the driven spherical p-spin [26, 28] are 



d t C(t,t') = -n(t)C(t,t') + J dt"T,{t,t")C(t",t') + J dt"D(t,t")R(t',t") , 

d t R(t, t') = -n(t)R(t, t') + J dt"T,(t, t")R{t", t') + 5(t - t') , (94) 

H(t)=T + J dt'[D(t,t')R(t,t') + T,(t,t')C(t,t')] , 



with the vertex and self-energy 

D = + a 2 ^- 1 =D + a 2 D, , 

^= l - P {p-l)RCf- 2 = RD' (C) , 
respectively. The dynamics can also be described with a "single-spin" Langevin equation of the form 



(95) 



<p(t) = -n(t)<p(t) + J dt' E(t, t'Mt') + p(t) , 

(p(t)p(t'))=2TS(t-t') + D(t,t>) . 

Note that S and D verify the detailed balance condition. From the expressions (95), one can rewrite Eq. (96), in the 
stationary case, in the following way 



<p(t) = -M*) + J dt' £(i - t')ip(t') + p{t) + ah(t) , 

(p(t)p(t')) = 2TS(t - t') + D Q (t - t') , (97) 
(h(t)h(t'))=D 1 (t-t') , 

where p(t) and h(t) are two uncorrelated Gaussian variables. Note that S and Do still depend implicitly on a as one 
has to solve the self-consistency equations (94) for C and R and substitute the result in S and D n . 

In the absence of drive and interpreting E and D as the response and correlation of a bath, its frequency-dependent 
effective temperature is 

T(u>) = V{UJ) = 2 1±^M {m 
K ' 2Re ff H 2Re[l + £(w)/(iw)] ' 1 ' 
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Given that S and D depend on R and C, if one finds that R and C are related by the fluctuation dissipation theorem, 
R(t) = —(3&(t)C(t), from the relation S = RD' (C) it follows that T(lu) = T and the bath is in equilibrium, as 
expected. If R. and C do not verify the fluctuation dissipation theorem, T(J) ^ T. 
Switching on the external drive one can compute its dissipated power 

W(t) = ah(t)ip(t) (99) 

and its average. One finds 

/>oo 
dt' R(t - t')[p(t') + ah(t')]) 
-oo 

f°° k f°° 

= (W)=a 2 dtR{t)D 1 {t) = a 2 - dt R^C^ 1 (t) , (100) 
Jo 2 J Q 

consistently with the result of [28] where the average of the injected power was explicitly computed for this model. 
One can prove that the entropy production <S* SS obtained from the rate 



diss 
t 



a I di!T-\t-l!)[h(t)<p(l!) + h{l!)<p{t)] , (101) 

J — oo 



verifies the fluctuation relation. Indeed, first rewriting the dissipative entropy production as 

n= — oc v ' 

and using the solution for ip n , in the generic notation of the previous section, one finds 

(103) 



where v, p and g are the Fourier transforms of the p — p correlator, the h — h correlator and the time-integrated X, 
respectively. It is easy to check that (j)diss{^) = 4>diss(l — A) and the fluctuation relation is then verified. 

Once again, we showed that the dissipative contribution to the entropy production satisfies the fluctuation relation 
with a modified temperature. Note however that this solvable example is non-trivial for at least two reasons. It is a 
clearly an-harmonic problem since X and D are themselves functions of C and R. The single spin ip is coupled to an 
equilibrated bath at temperature T and a self-generated "bath" at a different temperature. The temperature entering 
the fluctuation relation involves both, through the definition (98). This temperature is the one that one would observe 
by measuring the fluctuation-dissipation relation on the variable tp. 



C. The adiabatic approximation 



When a simple system is coupled to a complex bath with two (or more) time scales these are induced into the 
dynamics of the system. When the time-scales are well separated, an adiabatic treatment is possible in which one 
separates the dynamic variables in terms that evolve in different time-scales (dictated by the baths) and are otherwise 
approximately constant. 

In this subsection we use an adiabatic approach [10] to treat simple problems coupled to baths that evolve on 
different scales. The motivation for studying this type of problems is that they resemble glassy systems although in 
the latter the separation of time-scales is self-generated. 

We study the PDF of S^ tss and S T . The former satisfies the fluctuation theorem (at least in the harmonic case 
since corrections due to might appear in non-harmonic problems). We check that the adiabatic approximation 
does not spoil this feature. The latter, instead, does not satisfy the fluctuation theorem in general. However, it 
is interesting to understand under which conditions it satisfies the fluctuation theorem approximately. Indeed, in 
numerical simulations and experiments it has been customary to measure the dissipated power W t — ah{<p t ]<pt and 
then define the entropy production rate as a t — W t /T (where T is the temperature of the fast bath) . This corresponds 
to the definition of of , given in Eq. (85) for the harmonic oscillator. We show that when the particle is coupled to a 
bath that evolves in well-separated time-scales of does not satisfy the fluctuation theorem. 
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FIG. 5: Power spectrum of the entropy production rate (full line) as a function of the frequency for the harmonic oscillator 
with a = 0.5, k = 1, T/ = 0.6, T s = 2, 7/ = 1, 7 S = fcT s and r s = 250. The dot-dashed line is the "slow" contribution of H t , 
the dashed line is the "fast" contribution of wt- 



Let us consider again the Langevin equation (87) with V(|z| 2 ) = k\z\ 2 . In this case, the correlation functions can 
be calculated explicitly, see Appendix D. In Fig. 3 we report the autocorrelation functions, C(t) = (z t zo), for a = 0.5, 
k = 1, Tf = 0.6, T s = 2, 7/ = 1, 7 S = kr s and different values of r s . Clearly, for kr s = -f s ^> 7/ two very different 
time scales -related to the time scales of the two baths- are present. From the plot of Fig. 4 one sees that in the case 
kr s = 250 7/ the function x(C) is a broken line with slope —1/Tf at large C (short times) and — 1/T S for small C 
(large times). 

We want to show that, in this situation, the variable zt can be written as the sum of two quasi-independent 
contributions. Using the construction introduced in [10] we rewrite the equation of motion (87) as 



The variable h t is "slow" ; if we consider it as a constant in the first equation, the variable z t will fluctuate around 
the equilibrium position Zh — h/(k + "/ s /t s — ia) = H. The latter will -slowly- evolve according to the second 
equation in (104), in which we can approximate z t > ~ H±>. Defining the -fast- displacement of z t with respect to H t , 
Wt = z t — H t , we obtain the following equations for (wt, H t ): 



In this approximation, z t — H t + Wt is the sum of two contributions: w t is a "fast" variable which evolves according 
to a Langevin equation with the fast bath only and a renormalized harmonic constant fc + 7 s /r s , while H t is a "slow" 
variable which evolves according to a Langevin equation where the slow bath only appears. In both equations the 
driving force a is present, thus we expect both H t and w t to contribute to the dissipation. Note that w t and H t are 
completely uncorrelated in this approximation. 



1. The harmonic model coupled to two baths 




(104) 




(105) 
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FIG. 6: The function [Ce(p) - (e(-p)]/a+ for the harmonic oscillator with a = 0.5, k = 1, 6 = T) = 0.6, T s = 2, 7/ = 1, 
7 S = fcr s and r s = 1, 50, 250. The dashed line (FR) is a straight line with slope 1, the dotted line has slope Tf/T s = 0.3. 



2. The "potential" entropy production rate o~Y 

In the adiabatic approximation the term in equation (41) should become a boundary term. Indeed, the function 
T _1 (t), in the adiabatic approximation, becomes 



T-\t) = -6(t)+T-\t) , 



T 



(106) 



where the function T s 1 (i) is "slow", see e.g. Eq. (91). Inserting this expression in af > the first term gives a total 
derivative. The second term gives 



dt'T- l (t-t') 



dV dV 
^{t)^W)]+f a ^) — [r{t)] 



(107) 



Due to the convolution with the "slow" function T~ l {t), the fast components of r are irrelevant in the integral, while 
for the slow ones it is reasonable to replace r a {t) with f a (t') on the scale t s over which T^ 1 ^) decays. Thus one 
obtains a total derivative times the integral of T s _1 (f) which is a finite constant. Obviously this is not a rigorous proof 
and should be checked numerically in concrete cases. 



3. The "dissipative" entropy production rate a t tss 

The entropy production rate defined in Eqs. (80) and (89) can be rewritten in terms of H t and Wt- Recalling that 
T _1 (t) is defined by Eq. (91) one obtains (the details of the calculation are reported in Appendix E) 



diss 



2 aim 



w t w t 



H t H t 
T s 



(108) 



neglecting terms that vanish when af lss is integrated over time intervals of the order of t s . This is exactly the entropy 
production expected for two independent systems. 

To check that this approximation works well, let us introduce a "power spectrum" a{ui)dw as the contribution 
coming from frequencies [u,u + dcu] to the average entropy production rate, er* ss — duja(uj). From Eq. (84) we 
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get 



diss 

= 



d<t> 



diss 



dX 



2tt 



A=0 

2aujv{u)) 

TdRF 



dw 2awv(uj) 
2^ \D{uj)\ 2 



2auv(- 



\D(-u)\< 



OiU)v{u)) 



\D(-u;W I^HN 



(109) 



Substituting the expressions of v(u>) and of D(u>) appropriate for Eq. (87) we get the power spectrum o(u>) as a 
function of uj which is reported in Fig. 5 as a full line. The contributions of w t and H t , cr w (uj) and oh(u)), are 
obtained inserting in Eq. (109) the expression of v>(u>) and D(u) obtained from the two equations (105). They are 
reported as dashed and dot-dashed lines in Fig. 5. We can conclude that, for kr s 3> 7/, the adiabatic approximation 
holds and af lss — o™ + of , with o™ — 2almw t Wt/Tf and of = 2almH t H t /T s , and the two contributions are 
independent. Note that the average dissipation due to H is much larger than the one due to w. Finally, we can write: 



,(A) = <T(A) + (t> H {\) , <j> w ' H {X) = - lim t- 1 ln(cxp [ - Xo w ' H ]) . 



(110) 



Both 4> W (X) and 4> H (X) verify the fluctuation theorem, as the two equations of motion (105) are particular instances 
of the general case discussed in section III. The function Qi SS (p) is the Legendre transform of 4>diss(X) and will verify 
the fluctuation theorem. 



4- The "approximate" entropy production rate of 
In the same approximation, of is given, for 6 = Tf, by 

= ar + |V, and <j> Tf (A) = <T(A) + cj> H (XT s /T f ) ; (111) 

the contribution of H t is weighted with the "wrong" temperature, i.e. the temperature of the fast degrees of freedom. 
Indeed, as we have already discussed, (f>e(X) does not verify the fluctuation theorem. The function f(p) = [(e(p) — 
(e(—p)]/&+ [obtained from Eq. (86)] is reported in Fig. 6. As already discussed in section IV, when the time scales 
of the two baths are of the same order, kr s <~ 7/, the two baths act like a single bath at temperature 6 £ [Tf,T s ] and 
the function f(p) is approximately linear in p with slope X G [T//T s , 1]. When the time scales are well separated, 
kr s S> 7/, the adiabatic approximation holds; and one finds that f(p) has slope ~ 1 for small p and Tf/T s for large p 
(see Fig. 6). 

The results for kr s ~ 7/ are consistent with the ones in [41] where only the situation in which the two time scales 
are not well separated could be investigated. Indeed, when the dynamics becomes very slow the observation of large 
negative fluctuations of the entropy production requires a huge amount of computational time and the function f(p) 
can be calculated only in a narrow range of p around p = 0. Note that the value of p at which the slope of f(p) crosses 
over from 1 to Tf/T s depend on the values of the parameters a, 7/, 7 S , r s , etc. and can be of the order of 5, while in 
numerical simulations one can usually reach values of - at most - p <~ 3. Thus, the observation of curves like the one 
reported in Fig. 6 in numerical simulations of glassy systems is a very difficult task. 



2alm 



w t w t + H t H t 
Tf 



VI. GREEN-KUBO RELATIONS 



It was proven in [32, 33] that the fluctuation theorem implies, in the equilibrium limit (07 — > 0) the Green-Kubo 
relation for transport coefficients. This is a particular form of the fluctuation-dissipation theorem. In this Section 
we discuss how one links the modified fluctuation theorem - in which we replaced the external bath temperature by 
the (frequency dependent) effective temperature of the unperturbed system - to the modification of the fluctuation 
dissipation relation. 



A. General derivation 



Let us recall briefly how the Green-Kubo relation can be obtained from the fluctuation theorem. Suppose that 
we apply a (constant) driving force E to a system in equilibrium. This generates a corresponding flux J t (e.g. if E 
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is an electric field Jt is the electric current) such that, close to equilibrium, the dissipated power can be written as 
W t — EJ t . The entropy production rate is then 



are the moments of <S T , i.e. 



The fluctuation theorem can be written as 0(A) — 0(1 — A) where 0(A) is defined in Eq. (25). The derivatives of 0(A) 

= (-l)*" 1 hm t-^c , (113) 



_ d fe 



A=0 T ^°° 



where (A 2 ) c = (A 2 ) — (A) 2 and so on. Thus, 0^ ~ er+, and close to equilibrium (er + <~ 0) 0(A) is well approximated 
by a second order polynomial (corresponding to a Gaussian pdf), 

0(A)~0 o + 0iA+^A 2 . (114) 

Then the fluctuation relation, 0(A) = 0(1 — A), implies 2 = — 20i; from Eq. (113) and using time-translation 
invariance, 

02 = -201 => (7+ = / (cr t cr ) c • (115) 

Jo 

Substituting a t = EJ t /T one obtains 

E f 00 

(J)e dt (JtJo) E =o + o(E 2 ) , (116) 

that is to say, the Green-Kubo relation. 

Note that even out of equilibrium one can define a flux J t using a t as a "Lagrangian" , see Ref. [44] : 

Close to equilibrium a t is given by Eq. (112) and J t — Jt/T. If, in the absence of a drive, the system has a non trivial 
effective temperature, the entropy production rate should be defined as in Eqs. (80) and (89). Then the flux J t is 
given by 

Q diss rt rt 

Jt = -^=4lm J dt'T-\t-t')z t z t ' =2 J dt'T- l {t-t')[y t x tl - x t yv\ ■ (H8) 
The fluctuation theorem for a dlss implies then a Green-Kubo relation for J t : 

poo 

{J) a =a dt (J t Jo) a =o + o(a 2 ) . (119) 
Jo 

The physical meaning of the latter relation becomes clear if one writes the flux J t in the adiabatic approximation 
discussed in the previous section; from Eq. (108) 



Jt = 21m 



w t w t , H t H t 



jw jH 

%- + , (120) 

If Is 



and Eq. (119) becomes 

(J w )a , (J H ) a _ a 



Tf 



I Jtl \ I'OC p oo 

+ ~Tr = TfJ dt ( J ? J o)«=o + ^J o dt(J t H J H ) a=0 + o(a 2 ). (121) 



Indeed, in the adiabatic approximation the Green-Kubo relation holds separately for Jf (with temperature Tf) and for 
(with temperature T s ). Eq. (119) encodes the two contributions and holds even when the adiabatic approximation 
does not apply and the contributions of the "fast" and of the "slow" modes is not well separated. 
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Note that the "classical" Green-Kubo relation involves the total flux J t = Jf + J t . For the latter one has, in the 
adiabatic approximation, 

poo />oo 

(Jt)a = (J?)* + (J t H )a = 7fr dt (Jf J -) a=0 + - / dt (J t H J H ) a=0 

1 f JO 1 s Jo 



= a 



dt 



{Jt Jo)a=0 (J{ 1 JQ H )q=0 



Tf 



~ a 



(122) 



dt 



Teff(t) 



(JtJo)a=0 



The latter relation is the generalization of the Green-Kubo formula that comes from the generalized fluctuation- 
dissipation relation discussed in the Introduction and Section III. It is closely related, but not equivalent, to Eq. (119). 



B. The Green-Kubo relation for driven glassy systems 

Eqs. (119) and (122) cannot be applied straightforwardly to driven glassy systems as for these systems the correlation 
function (JtJo)a is not stationary at a — at low temperatures. Indeed, the relaxation time of the latter grows very 
fast as a — > and at some point falls outside the experimentally accessible range: the system will not be able to reach 
stationarity on the experimental time scales and will start to age indefinitely. 

However, let us consider again the equation of motion (10) for a^O, where we assume that a stationary state is 
reached, 

7<i>t = -M« Vt+ dt' E a (i - t') <p t > + Pt + ah[ip t ] , (ptPr) = D a (t - t') . (123) 

J — oc 

The functions E a (f — t') and D a (t — t') depend (strongly) on a; indeed, as the term explicitly proportional to a is a 
small perturbation for a ~ 0, the main contribution to the a-dependence of the dynamics of tpt will come from the 
a-dependence of S Q and D a . If we compare the latter equation with Eq. (87), we see that setting a = in Eq. (87) 
is equivalent to setting a = without changing the functions E and D in Eq. (123). This will not affect too much the 
correlation function (JtJo)a if ot is small. Finally, we can write, for small a, 

r°° i 

(Jt) a ~a dt——U t J ) a , (124) 

^0 Teff(t) 

even if the limit a — > is not well defined. An analogous relation will be obtained from Eq. (119) (which is equivalent 
to the fluctuation theorem in the Gaussian approximation) within the same approximation. The latter relations can 
be tested in numerical simulations as well as in experiments. 



VII. SLOW PERIODIC DRIVE AND EFFECTIVE TEMPERATURE 



In this Section we discuss a means to measure an effective temperature associated to slow timescales of a non- 
equilibrium system by using the modification of the fluctuation theorem. 

A lesson we learn from the previous calculations (see e.g. Fig. 5) is that the work done at large frequencies is 
overwhelmingly larger than that done at very low frequencies - precisely the one we wish to observe in order to detect 
effective temperatures. One way out of this is to choose a perturbation that does little work at high frequencies: a 
periodically time-dependent force that derives from a potential cos(£lt)V (r) , with 1/fl of the order of timescale of the 
slow bath t s . Let us show this for a one dimensional system, the generalization being straightforward. 



A. General derivation 



Let us consider a single degree of freedom r moving in a time- independent potential V(r) and subject to a periodically 
time-dependent field cos(Qt)V(r), and in contact with a 'fast' and a 'slow' bath with friction kernel, thermal noise 
and temperature (p?,gf,Tf) and (p s ,g S7 T s ), respectively: 

dt' g f (t-t')r(t')+ J dt' g s (t-t')r(t') = -^[r(t)]+pf(t)+p^t)-cos(nt)^[r(t)} , (125) 
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The time scale of the time dependent field 1/0 is of the same order as the one of the 'slow' bath. The work in an 
interval of time (0, r) done by the time-dependent potential is: 

r dv r 

W T = - cos(fii') — rdt' = -V(t) + V(0) + n / sin(Qi') V dt' . (126) 
Jo c* r Jo 

Only the last term grows with the number of cycles, so for long times we can neglect the first two. Now, integrating 
(125) by parts, we obtain 

mr(t) = -f_ dt' g f (t-t')f(t')-^[r(t)}+pf + h(t) - h(t)^[r(t)} (127) 

h(t) = - [ dt' g s {t-t')r(t') + p s (t) . (128) 

J — oo 

where h(t) — cos(fii). In the adiabatic limit when both the timcscalcs of the slow bath and the period l/fi of the 
potential V are large, h(t) and h(t) are quasi-static. Hence, r has a fast evolution given by Eq. (127) with h and h 
fixed and it reaches a distribution [10] 

-P J (v+hV+g J (0)4~hr) 

P(r/h, h) = - - r . (129) 

J dr e -0f(v+hV+ gf (O)4-hr) 

The denominator defines Z(h, h) and F(h, h) = —fij 1 \nZ(h, h). Note that F(h, h{t)) is periodically time-dependent 
through h. The approximate evolution of h is now given by Eq. (128) with the replacement of r in the friction term 
by its average dF gl' h ^ with respect to the fast evolution: 

h{t) = J t dt' gs (t-t') dF ^ h h \ t') + P s (t). (130) 

Eq. (130) is in fact a generalized Langevin equation for a system coupled to a (slow) bath of temperature T s . Indeed, 
it can be shown [10] to be equivalent to a set of degrees of freedom yi evolving according to the ordinary Langevin 
equation: 



d d o 

m ^ + ^dt +n > 



OF 



with (^i{t)^j(t')) — 2T s 'jj5ijS(t — t'), provided that the Fourier transforms g s (u) an v s {lo) of friction kernel and noise 
autocorrelation can be written as: 

ffM = E— 7 % =T ^)=2T S ^— — ljA \ (132) 

where are the roots of —uijuJ 1 + i^jUj + Qj = 0. 

Within the same approximation leading to (130), the average of V(r) over a time window A that is long compared 
to the short timescale, but sufficiently slow that we can consider that h and h are constant is 

( t+A V(r(t')) dt' <~ A / dr P(r/h, h) V(r) = A ^^' ^ (133) 
Jt J dh 

so that we obtain 

w^nf M nt') dF{hMt/)) dt' = - f T d J^!pi d , (134 ) 

Jo dh Jo ot 

which tells us that for long time intervals the work done by the original time-dependent potential V is indeed the 
same as the work done by the time-dependent effective potential F in (130). 

The fluctuation theorem then holds for the distribution of this work, with a single temperature T s . We conclude that 
the distribution of work due to a slow perturbation satisfies the fluctuation theorem with only the slow temperature, 
and can be hence used experimentally to detect it. 
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B. Experimental realization 

The simplest application of the above general result is obtained considering V(r) = hr and V(r) = kr 2 . Then, 
grouping together the two noises in a single one with friction g = gf+g s and correlation v = TfVf + T s v s as described 
in Sect. Ill, Eq. (125) simply becomes 

mr(t) + [ dt' g(t - t')r(t') = -kr(t) + p(t) + h cos(Ot) . (135) 

J — oo 

This equation describes for instance the motion of a Brownian particle moving in an out of equilibrium environment 
and trapped by an harmonic potential whose center oscillates at frequency 0. A concrete experimental realization 
of this setting has been already considered in [35]: Silica beads of ~ 2/im diameter were dispersed in a solution of 
Laponite (a particular clay of ~ 30nm diameter) and water. The Laponite suspension form a glass for large enough 
concentration of clay and provides the nonequilibrium environment. The Silica beads are Brownian particles diffusing 
in such environment. They can be trapped by optical tweezing, and the center of the trap can oscillate with respect to 
the sample if the latter is oscillated through a piezoelectric stage. In [35] the mobility and diffusion of tracer particles 
were measured obtaining an estimate of T e ff(Q). Here we propose to measure the work done by the trap on the 
tracers. Indeed, the work dissipated in (0, r) is linear in r(t) so it should be possible to measure it simply through 
the measurement of r(t): 

W T = n~h[ dt' sm(nt')r{t') ; (136) 
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note that, as W T is linear in r(t), it is a Gaussian variable. With a simple calculation one finds 

r^oo 2{W T ) 2Reg{Q) tn ' V ' 

This means that the (Gaussian) PDF of S^ tss = W T /T e ff(il) satisfies the fluctuation relation. If the two baths are 
modeled as in Sect. IV with k-f s = t s 3> 7/ one has T e ff(Sl) = T s for f2r s < 1 [see Eq. (90)]. The measurement of the 
distribution of the work (136) allows for the measurement of T s . Note that other experimental settings described by 
the same equations should exist. 



VIII. CONCLUSIONS 



We studied the extension of the fluctuation relation to open stochastic systems that are not able to equilibrate with 
their environments. 

We used the simplest example at hand to test several generalized fluctuation formulas: a Brownian particle in a 
confining potential coupled to non-trivial external baths with different time-scales and temperatures. Independently 
of the form of the potential energies, due to the coupling to the complex environment, the particle is not able to 
equilibrate. Its relaxational dynamics is characterized by an effective temperature, defined via the modification of 
the fluctuation-dissipation relation between spontaneous and induced fluctuations [9]. When no separation of time- 
scales can be identified in the bath, the effective temperature is a non-trivial function of the two times involved. 
Instead, when the bath evolves in different time-scales each characterized by a value of a temperature, the two-time 
dependent effective temperature is a piece- wise function that actually takes only these values, each one characterizing 
the dynamics of the particle in a regime of times. 

Several authors discussed the possibility of introducing the effective temperature in the fluctuation theorem to 
extend its domain of applicability to glassy models driven by external forces [12-15]. After summarizing our results 
we shall discuss how they compare to the proposals and findings in these papers. 



A. Summary of results 



We here examined carefully different definitions of entropy production rate that are not equivalent when the effective 
temperature is not trivially equal to the ambient temperature. We found that: 
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1 . The PDF of the "dissipative" entropy production S^ tss that involves the frequency dependent temperature, 

Sf- = ^ £ — "^"^ K) = r dt f 72 dt>T-\t-t<)r a (t)KW)\ , d38) 

with T _1 (t) the Fourier transform of 1/T(u>), the effective temperature of the relaxing system, see Eq. (18), 
verifies the fluctuation theorem exactly for harmonic system, ft also holds for general systems connected to 
baths with different tcmperaures acting on widely separated scales. 

2. For nonlinear systems in contact to noncquilibrium baths acting on overlapping timcscalcs an additional term 
S^ , see Eq. (41), has to be included in the entropy production to make the fluctuation relation hold strictly. 
Our numerical simulations suggest that, surprisingly enough, the effect of this extra 'internal' term is even then 
quite small. 

3. The PDF of S® with af = W t /Q, W t the power dissipated by the external force and 6 a free parameter with 
the dimensions of a temperature, does not satisfy the fluctuation theorem in general for any choice of 9. 

The large deviation function, (p) , still shows some interesting features revealing the existence of an effective 
temperature. When the bath has, say, two components acting on different time scales and with different tem- 
peratures, the function [(&(p) — Ce( — P)]/&+ ma y have different slopes corresponding to these two temperatures, 
one at small p and the other at large p. The separation of time-scales of the bath translates into a separation 
of scales in the function [Ce(p) — Ce( — P)]/°+- 

When the time scales of the baths are not separated, and one records the large deviation function for not too 
large values of p only, the fluctuation theorem is verified approximately if O is suitably chosen. Note that the 
value of 9 found in this way is not equal to the effective temperature T e f / that enters the modified fluctuation- 
dissipation relation. Instead, when the time-scales are well separated, the two scales in the large deviation 
function arc clearly visible and a single fitting parameter is not sufficient to make the fluctuation theorem hold. 

4. If two time scales are present in the dynamics of a system and the applied perturbation is periodic with frequency 
ft < 1/t s , t s being the largest relaxation time, the PDF of the power dissipated over a (large) number of cycles 
verifies the fluctuation relation with temperature T s = T e //(fi). This is probably the easiest way of detecting 
the effective temperature by means of the fluctuation relation. 

These results should apply to driven glassy systems as discussed in section V and are indeed consistent with recent 
numerical simulations [41]. Models like the one discussed here have been recently investigated [10, 35-37] to describe 
the dynamics of Brownian particles in complex media such as glasses, granular matter, etc. Brownian particles are 
often used as probes in order to study the properties of the medium (e.g. in Dynamic Light Scattering or Diffusing 
Wave Spectroscopy experiments). Moreover, confining potentials for Brownian particles can be generated using laser 
beams [39] and experiments on the fluctuations of the power dissipated in such systems are currently being performed 
[35,47]. 

B. Temperatures 

It is important to summarize the different definitions of effective temperature we considered and the relations 
between them. We defined the effective temperature in the frequency domain in equation (18) as a property of the 
bath which can also be measured from the ratio between correlation and response functions in the frequency domain. 
As we discussed above, the same effective temperature enters the correct definition of entropy production rate in the 
frequency domain, see equation (138). Thus, experiments working in the frequency domain should observe the same 
effective temperature from the fluctuation-dissipation relation and the fluctuation relation. 

In the time domain the situation is slightly more complicated. On the one hand, the effective temperature obtained 
from the fluctuation-dissipation relation in the time domain, defined for example by Eq. (9), is not the Fourier 
transform of T(lu). A convolution with the correlation function is involved in the relation between T(w) and T e ff(t). 
On the other hand, the effective temperature T _1 (i) entering the entropy production is exactly the Fourier transform 
of 1/T(w), see again Eq. (138). This can give rise to ambiguities when working in the time domain. 

Most of these ambiguities disappear as long as the time scales in the problem arc well separated. In this case, on 
each time scale a well defined effective temperature can be identified, and this temperature enters both the fluctuation- 
dissipation relation and the fluctuation relation: see e.g. the curve for r = 250 in Fig. 4 and the expression of <S* SS 
in the adiabatic approximation, eq. (108). This is essentially related to the validity of the adiabatic approximation 
discussed in section VC. 



30 



The difference is relevant when the time scales of the two baths are not well separated, and a single effective 
temperature cannot be identified, see the curve for t = 1 in Fig. 4. In this case, we found that the fluctuation 
relation holds with -approximately- a single fitting parameter but this temperature is not clearly related to the 
fluctuation-dissipation temperature in the time domain. This is indeed what is observed in numerical simulations on 
Lennard-Jones systems [41]. 

Let us remark again that, when applying these results to real glassy systems in finite dimension, one should take 
care of the possibility that the effective temperature has some space fluctuations due to the heterogeneity of the 
dynamics [38]. The extension of our results to such a situation is left for future work. 

C. Discussion 

Several proposals to introduce the effective temperature into extensions of the fluctuation theorem have recently 
appeared in the literature. Let us confront them here to our results. 

Sellitto studied the fluctuations of entropy production in a driven lattice gas with reversible kinetic constraints [12]. 
When coupling this system to an external particle reservoir with chemical potential ji, a dynamic crossover from a 
fluid to a glassy phase is found around n c . The glassy nonequilibrium phase is characterized by a violation of the 
fluctuation dissipation theorem in which the parametric relation between global integrated response and displacement 
yields a line with slope fj, e ff [48]. 

One drives this (possibly already out of equilibrium) system by coupling two adjacent layers of the three dimensional 
periodic cube to particle reservoirs at different chemical potentials, fi + and The former is allowed to assume value 
corresponding to the glassy phase, n+ > fi c , while fi- is always below \i c . The results of the Montecarlo simulation 
are consistent with a generalized form of the fluctuation-theorem: 

a T = Jri/J-eff ~ M-) i ( 139 ) 

where o> is the entropy production, J T is the particle current in the direction of the externally imposed chemical 
potential gradient averaged over a time- interval of duration r; l^eff is an effective chemical potential and is the 
chemical potential of one of the layers. When the chemical potentials of the two reservoirs are in the fluid phase, 
fi e ff = /i+ and the usual fluctuation relation holds. Instead, when fi + is in the glassy phase, Sellitto found that 
Eq. (139) holds with /i e // taking the value appearing in the violation of fluctuation-dissipation theorem in the aging 
regime of the undriven glassy system at /i+ . 

The formula (139) differs from the ones we found to describe the oscillator problem in that in our case, when 
translating from temperature to chemical potential, the full time-dependent fi(t) enters. Strictly, we believe that 
this improved definition should also apply to the lattice gas model. However in the case studied by Sellitto the fast 
dynamics is an "intra-cage" dynamics that likely does not contributes to the current. This is a case in which the 
perturbation does not produce dissipation at high frequency so that the difference arising from fi(t) ^ fj, e ff should be 
tiny in this case (see Sect. VII). 

More recently, Crisanti and Ritort [14] found that the probability distribution function of the fluctuations of heat 
exchanges, Q, between an aging random orthogonal model in its 'activated regime' (a long-time regime in which the 
energy-density decays as a logarithm of time) and the heat bath is rather well described by a stationary Gaussian part 
and a waiting-time dependent exponential tail towards small values of Q. Assuming that these events are of two types 
('stimulated' and 'spontaneous') they proposed to fit the ratio between the PDF of positive and negative 'spontaneous' 
Q's in the form of a fluctuation theorem, i.e. to be proportional to e~ 2( 2/ A , and relate A to the effective temperature of 
the fluctuation-dissipation relation. They found good agreement. Crisanti, Ritort and Picco are currently performing 
simulations to test this hypothesis in Lennard-Jones mixtures [49]. 

Another development is an attempt to generalize the situation considered by Crooks. He considered a problem 
that starts from equilibrium in zero field and evolves according to some stochastic dynamic rule in the presence of 
an arbitrary applied field [50] and found that the ratio between the probability of a trajectory and its time-reversed 
one is given by e~Pf° max dth(t)0(t) w j^ n /j^) ^ ne time-dependent external field that couples linearly to the observable 
O. For simplicity, let us focus on O = (f> with (f> a scalar field characterizing the system. In [13] the extension of this 
relation to the initial non-equilibrium 'glassy' case was conjectured. Separating the external fields h and <j> in their 
fast and slow components [51], h = hf + h s and cj> — tfif + <p s , one then proposes that the PDFs of the trajectories of 
the slow components satisfy a relation similar to Crooks' with the temperature replaced by the effective temperature 
(for a glassy non-equilibrium system with two correlation scales [3]). This is indeed very similar to what we have done 
in this paper. 

Finally, let us mention the work of Sasa [15] where he introduces an effective temperature in his definition of entropy 
production for the Kuramoto-Sivashinsky equation. 
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APPENDIX A: DIRICHLET BOUNDARY CONDITIONS FOR THE WHITE BATH 



A second possibility to calculate the functional integral in Eq. (68) is to impose Dirichlct boundary conditions 
z(— t/2) = 2(t/2) = 0. However, in this case it is possible to calculate <j)(X) only for m = 0. The distribution of z t is 
obtained substituting p u = D{ld)z u in Eq. (72): 



1 ( x du 2 _ 



exp 



1 o 1 c\ ^ 2^1 

dt zt [ k +a - %ia~i— - 7 ^2 ) z * 



V[z t ]= exp 
From Eq. (68) 

/If 00 / d d 2 \ 

dz t exp -^fj dt z t ik 2 + a 2 -2ia 1 [l-2\xr(t)]--l 2 ^U t 



(Al) 



(A2) 



where Xr(t) is the characteristic function of t £ [—t/2, t/2]. At the leading order in t, as the correlation function of 
Zt decays exponentially on a time scale To = 7& -1 , we can integrate out the portion of the trajectory that is outside 
the interval [—t/2, t/2] both in the numerator and the denominator, and we simply obtain 



f [ 1 C r/2 ( d d 2 \ 

(exp[-AS T ]) = AT 1 J dz t exp -— J dt z t [k 2 + a 2 - 2w 7 (l - 2A)— - 7^ J «t 



(A3) 



We have then to find the eigenvalues of the operator appearing in the integral. This corresponds to find the solution 
of the equation 



Jz t 



d d 2 \ 

k 2 + a 2 - 2i ai {l -2\)-- 1 2 —\z t = Ez t 



(A4) 



with boundary conditions z(t/2) — z(—t/2) = 0. Note that the operator J is Hermitian, thus the eigenvalues are 
real; they are given by the following expression: 



E n (X) =fc 2 +4a 2 A(l-A)+ r 



2^2 



(A5) 



with n = 0, 1, • • • . For each n the integration is performed on one complex variable and we get 



r 1 rr/2 

(exp[-A5 T ]) =M 1 / dz t exp \ dt z t Jz 

J ^l 1 J -t/2 



E n (0) 
■ \ E n (X) 



n 



(A6) 



n=0 



recalling that the constant Af is simply the numerator calculated in A = 0. Finally we obtain, defining w — mr/r, 



<t>(\) = lim - V In 

T — >OD T ' J 



E n (\) 



-oct— E n (0) 



In 



4a 2 A(l - A) 



k 2 + ~f 2 UJ 2 



(A7) 



The latter expression verifies obviously the fluctuation theorem. Moreover, in the m — case Eq. (74) is equal to 
Eq. (A7), as one can check using suitable changes of variable in the integral. In this simple case, £(p) can be computed 
exactly. Starting from Eq. (A7) one has 



0'(A) 



dio 



4a 2 (l - 2\) 



2a 2 (l-2A) 



o 7r 7 2 w 2 + k 2 + 4a 2 A(l - A) jy/k 2 + 4a 2 A(l — A) ' 



(A8) 
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and, recalling that (f)(0) = 0, 



0(A) = / dfi <p'(fi) = j- 1 [^k 2 + 4a 2 A(l - A) - k] . 



The function £(p) is defined by 



C(p) = mm[\pa + - 0(A)] - A* p( r + - 0(A*) , 



(A9) 



(A10) 



where <r + = 2a 2 /(jk) and A* is defined by 0'(A*) — pcr+; hence, 



fc(l - 2A*) 



V / fc 2 + 4a 2 A*(l- A*) 



1 -p* 



a 2 p 2 + fc 2 



and finally 



C(p)= 7 - 1 U- 



a 2 p 



1 - m 



/ a 2 + fc 2 
a 2 p 2 + A: 2 



/ a 2 + fc 2 
a 2 p 2 + fc 2 



From the latter expression it is easy to verify that 



C(p) - C(-p) = 



2a 2 p 
k 



(All) 



(A12) 



(A13) 



as stated by the FT. Defining r = 7/fc, the relaxation time of the correlation function of z t , and <r — <7+t /2 = a 2 /k 2 , 
the (adimensional) entropy production over a time to/2, we obtain 



C(p) = V 1 [l +p<T - V / (1 + ^0)(1+P 2 ^0) 



(A14) 



APPENDIX B: FLUCTUATION THEOREM FOR MANY EQUILIBRIUM BATHS AT DIFFERENT 

TEMPERATURES 

We compute the function 0(A) in the case in which the driven oscillator is coupled to N colored baths with generic 
memory functions and in equilibrium at different temperatures. The violation of the fluctuation-dissipation theorem 
for the relaxing particle in such an environment was discussed in Sect. Ill A 2. As discussed there, the equations 
are mathematically equivalent to the ones discussed in Sect. Ill; thus the strategy as well as many details of the 
calculation are the same as in this Section. 

The equations of motion are 



N -oo N 

+ ds gi(t- s)z s = -KZt + ^pn , (Bl) 

i=i J -°° »=i 

with k = k — ia. The thermal noises satisfy 

(pitpjo) = (PitPjo) = , 
(Pitpjo) = 5 lj T i v i {t) . 

By causality, the functions g%[t) must vanish for t < 0. As the baths are in equilibrium at temperature T i: the 
functions Vi(t) and #,(£) are related by Eq. (55): 

Vi {t) = Tilgiit) + 9i (-t)] = T i9i {\t\) , fB3 , 
T i9i {t) - 9(t)ui(t) . 1 ' 

In the frequency domain Eq. (Bl) becomes 

_ Si _ Si 

w —mu) 2 + n — iuLi^ i gi(uLi) D(u>) ' 
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where D(u>) = —mu 2 + k — iw gi{u>). 



The dissipated power is given by 

/oo _ 
ds gi(t - s)z t z s + 2Re ^ z t p it = W t - ^ W it 
-oo 



dH 
~dT 



(B5) 



where as in the previous cases Wt — 2almz t z t is the power injected by the external force and Wt — 2Re ds gi(t — 

s)ztz s — 2Keztf>it is the power extracted by the i-th bath. 

The first definition of entropy production rate, Eq. (85), gives (in the following, we will always substitute 
_jy^oo , r°° dui 

2ff Z^n= — oo J — oo 2ff 



Aw ' S^L-oo — * J^oo Iff as tne error is C(l) for t — > oo, see section III): 



rfcj 2c 



2tt e 



(B6) 



Substituting — Pito I D{uj), we obtain 

//■OO j 
^ cx P -y ^X>*4M* 



(B7) 



where A x (uj) is a, N x N real matrix which elements are given by 



4M 



A 2aw 



Tii/i(w) |£>H| 2 O ' 



Then, 



CH 7 ™ 11 dot A° K) ^ 27T 



In 



det A x (lu) 



det A°(w) 



(B8) 



(B9) 



The determinant of a matrix of the form A A = c f (Jy + A& satisfies the relation 

det A A 1 w ^ 
d^ = l + A6^ Ci ; 

we finally obtain 



o ( _ )( A)= / g In 



1 - 



(BIO) 



(Bll) 



In general, it does not exist a choice of Q such that 4>e(\) verifies the fluctuation theorem, i.e. <pe(X) ^ <^e(l — A). 

For the second definition, given by Eq. (80), the computation is identical to the one of the previous section with 
the substitution — > T(u>), where T(u>) is given by Eq. (61). The result is then 



fe s (A) = r ^ m fi - 2aXuJ ^ T Tt ] ] = f - m 

^"^ ' 27T [ Te//M|#M| 2 J 2lT 



(B12) 



Observing that 



D(-w) = D(lu) - 2ia , 

\D(-lu)\ 2 = \D(lu)\ 2 - 2aiuJ2Mu) 



(B13) 



and it is now easy to show that 4>diss(X) — fidissO- — A). 
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APPENDIX C: ENTROPY PRODUCTION OF THE THERMAL BATHS 



We will discuss here a different definition of entropy production rate based on the power extracted by the thermal 
bath instead of the one injected by the driving force. The two differ by a total derivative if there is only one bath, so 
their asymptotic distributions should be identical if boundary terms can be neglected. However Van Zon and Cohen 
[42] showed in a particular case that this argument is not correct, see also [43]. 

If there are many baths equilibrated at different temperature, the study of the entropy production extracted by 
each bath allows to separate the different contributions to the total entropy production weighting each one with the 
right temperature. We will first discuss the case of a single bath, and later we perform the computation for the general 
case. We will show that the entropy production rate defined in this way satisfies the fluctuation theorem as outlined 
in [44] ; unfortunately, in our computation we neglect all boundary terms so we cannot check if these terms modify the 
asymptotic distribution, as observed in [42, 43]. Anyway, the modification can be proven to be eventually relevant 
only for \p\ > 1 [42, 43], so the results we will discuss should hold at least for \p\ < 1. 

1. One bath 

The entropy production rate of the bath is defined as: 

a bath = t = - 2Rez t p t ] = at+ 0*?L, (CI) 

where Wt has been defined in section III and is the power extracted from the system by the thermostat. Since a t and 
a t differ only by a total derivative [see Eq. (CI)], for large r we have 

S bath = St + 0[h( t /2) - H{-t/2)\ =S t + /3AH . (C2) 

The first term in the r.h.s. is 0(t) and has fluctuations 0(y/r), while the second term has zero average and its 
fluctuations are also 0(1). On this ground one usually neglects the second term and concludes that 

C b ath(p) = C(p) ; (C3) 

the two definitions of entropy production rate arc equivalent and both distributions verify the fluctuation theorem. 
Van Zon and Cohen studied a model very similar to the one considered here, but where the distribution of S T is 
Gaussian while the distribution of AH shows exponential tails. In such a case, the large fluctuations of S^ ath are 
dominated by the distribution of AH even if ((S T - Tcr+) 2 ) > ((AH) 2 }, see Ref. [42, 43] and references therein. The 
modification of (bath(p) occurs for \p\ > 1; then, for |p| < 1, Cbath(p) still verifies the fluctuation relation. 



2. Many baths 



If the system is coupled to many equilibrated baths, in addition to the definitions of entropy production rate given 
by Eqs. (85) and (80), a generalization of Eq. (C2) is possible, given by the sum of the power dissipated by each bath 
divided by the corresponding temperature: 



baths 



N TJ7 



(C4) 



This quantity takes into account heat exchanges between the baths, and its average value does not vanish at a = 0, 
as we shall see in the following. 

Let us compute <f>bathsW- We have from Eq. (C4): 
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2 , 
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^3 Tj 
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ij 



IUJ 



D(uj)T 3 D(uu)Ti 



(C5) 
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If wc define the functions 



p(uj) = iu;D(uj) , 



(C6) 



we obtain 



/pOO J 
dPi* exp - J ^Y,P^ A U")Pi 



where A x (lo) is a N x N matrix which elements are given by 



(C7) 



X 



T iVi {uj) \D{u)\< 



F(uj) + 



p(u) p(uj) 



(C8) 



Then, 

(f>baths (A) 

The matrix A has the following form: 
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T^OO 
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(C9) 



(CIO) 



where (j, = X/\D(uj)\ 2 , c t = T^lj) and hj = F(uj) + ^ 
the following form: 



Its determinant is an order N polynomial in /i of 



dct A x 
dct A° 



i<j<k 



bii bij 



bik 
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Oji Ojj o jk 
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(Cll) 



Let us compute the coefficients explicitly. We will define Ty by TV 1 = T^ 1 — T- 1 . The coefficient of A 2 is given by 
a sum of determinants of the form 



F + T, 
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Ti 
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Ti Tji 



2 ' 



(C12) 



where we first subtracted the first column to the second column, and then subtracted the second row to the first row. 
We want now to show that all the coefficients of the higher powers of A vanish. Consider for example the coefficient 
of A 3 . It has the form 



-F I rp I rp F I rp I rp \ rp I rp 

* + Tk T~ F + + —. ^ + + jr 
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Ti 
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c 1 J P l JL JL JL. 
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(C13) 



where we subtracted the first column to the second and third column, and then subtracted the third row to the first 
and second row. The same argument applies to all the other coefficients up to order N . Finally, we get 
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and 



fibathsW 
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2^ 
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(C15) 



The first term in the logarithm is proportional to a and is related to the power injected by the external force, while 
the second term accounts for heat exchanges between the baths and does not vanish at a = 0. Finally, observing that 



D(-lo) = D{uj) - 2ia , 

\D(-lu)\ 2 = \D(oj)\ 2 -2aoj^Vi{w) 



(C16) 



and using the same trick we already used above it is easy to show that <j>baths{^) = 4>baths(l — A). Thus Cbathsip) 
should verify the fluctuation relation at least for \p\ < 1, if the contribution of boundary terms is not negligible. This 
result is of interest for the study of heat conduction and is similar to the one discussed in Ref. [44] . 



APPENDIX D: CORRELATION FUNCTIONS OF THE HARMONIC OSCILLATOR COUPLED TO 

TWO BATHS 



In the harmonic case, V(|z| 2 ) = jM 2 ' the correlation function of a variable z t , whose time evolution is given by 



Eq. (87), can be computed analytically [10]. In the frequency domain, Eq. (87) reads: 



pL + pI 
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where D{u) = k - iujjf - jj^r ■ Recalling that (p£p£,) = 47T7 / T / <5(w + J) and (p>*,) = jf^5(uj + u'), and 
defining C(u) from {z u z'J) = 2tt8{lo + uj')C(uj) we get 
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The function (1 — lot s )D{lo) is a polynomial in uo and its zeros are given by u — — i"f± where 
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2jfT s 
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and Re7± > 0. The response function is then given by 
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and the correlation function is given by 
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In the case a = 0, and in the limit 7/ « ] s « kr s where the time scales of the two baths are well separated, one 
obtains 
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and 



C(t) 



7~s 



i2(t) = 



Is) 2 

Is 



-t/r 



T f T s 

kr s + 7/ 
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If 



(D7) 



From the latter expressions it is easy to check that one has R(t) ~ —(3fO(t)C(t) for short times (t <C r s ) and 
i?(i) ~ —(3 s 6(t)C(t) for large times (t ~ r s ). The same behavior is found in the limit of small dissipation (small a), 
as one can check plotting the exact expression for the functions R(t) and C(t). 



APPENDIX E: THE EXPRESSION OF a s IN THE ADIABATIC APPROXIMATION 

We start from the expression (89) for af ss and from Eq. (91). Then (remember that ds S(t - s) = 6(0) = 1/2 
in our convention): 



ztz t , Is ( - T t x • 



dt'T-\t-t')z t z t , = ^+ - !L—— [i-±L) dt' e-W-^ZtZf . (El) 
2T f 2Q7> 7/ (t s ) 2 \ TfJ' 



We substitute Zt — Ht + wt and neglect all the terms proportional to HtWt'. indeed, such terms vanish when <j* ss is 

integrated over time intervals of the order of t s , as, on such time scales, (w t ) — while H is constant. The first term 
gives then 

ZtZt = H t H t + w t w t , E2 ^ 

2T f 2T f [ ' 

In the second term, as O ~ l/r s , we approximate dt' e -0 ^^* ^z t i ~ H t /Q to obtain 

^-^W t ~^fl-^UiJ t . (E3) 



2^ 2 T /7/ (r s ) 2 V T f J 2T S \ T f 

The (imaginary part of the) sum of these two terms times 4a gives Eq. (108). 
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